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Abstract. Extending investigations of Yarahmadian and Zumbrun in 
the strictly parabolic case, we study time-asymptotic stability of arbi- 
trary (possibly large) amplitude noncharacteristic boundary layers of a 

QQ , class of hyperbolic-parabolic systems including the Navier-Stokes equa- 

tions of compressible gas- and magnetohydrodynamics, establishing that 
linear and nonlinear stability are both equivalent to an Evans function, 
or generalized spectral stability, condition. The latter is readily check- 
able numerically, and analytically verifiable in certain favorable cases; in 
(-H I particular, it has been shown by Costanzino, Humpherys, Nguyen, and 

Zumbrun to hold for sufficiently large-amplitude layers for isentropic 
ideal gas dynamics, with general adiabiatic index 7 > 1. Together with 
these previous results, our results thus give nonlinear stability of large- 
amplitude isentropic boundary layers, the first such result for compres- 
sive ( "shock- type" ) layers in other than the nearly-constant case. The 
analysis, as in the strictly parabolic case, proceeds by derivation of de- 
tailed pointwise Green function bounds, with substantial new technical 

^t^ ' difficulties associated with the more singular, hyperbolic behavior in the 

fT^ . high-frequency/short time regime. 
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1. Introduction 

In this paper, we study the stability of boundary layers assuming that 
the boundary layer solution is noncharacteristic, which means, roughly, that 
signals are transmitted into or out of but not along the boundary. In the 
context of gas dynamics or magnetohydro dynamics (MHD), this corresponds 
to the situation of a porous boundary with prescribed inflow or outflow 
conditions accomplished by suction or blowing, a scenario that has been 
suggested as a means to reduce drag along an airfoil by stabilizing laminar 
flow; see Example 11.11 below. 

We consider a boundary layer, or stationary solution, 

(1) U = U{x), lim U{z) = U+, U{0) = Uo 

of a system of conservation laws on the quarter-plane 

(2) Ut + F{U), = {B{U)a,),, x,t>0, 

C/,F G M", S G M"^", with initial data U{x,0) = Uo{x) and Dirichlet type 
boundary conditions specified in ([5]), ([H) below. A fundamental question 
connected to the physical motivations from aerodynamics is whether or not 
such boundary layer solutions are stable in the sense of PDE, i.e., whether 
or not a sufficiently small perturbation of U remains close to U, or converges 
time-asymptotically to U, under the evolution of ([2]). That is the question 
we address here. 

Our main result, in the general spirit of [ZHl IMaZ3[ IMaZ4[ [Z3l iHZllYZ] . 
is to reduce the questions of linear and nonlinear stability to verification of 
a simple and numerically well-posed Evans function, or generalized spectral 
stability, condition, which can then be checked either numerically or by the 
variety of methods available for study of eigenvalue ODE; see, for example, 
[BnilB?2llB7ZllBDGl[lHZ2l[PZllFSllBim,Z.HLZi|HLyZlt|HLyZ2t[CH^^^ 



Together with the results of [CHNZj . this yields in particular nonlinear sta- 
bility of sufficiently large-amplitude boundary-layers of the compressible 
Navier-Stokes equations of isentropic ideal gas dynamics, with adiabatic 
index 7 > 1, the first such result for a large compressive, or "shock- type" , 
boundary layers. The main new difficulty beyond the strictly parabolic 
case of [YZ] is to treat the more singular, hyperbolic behavior in the high- 
frequency regime, both in obtaining pointwise Green function bounds, and 
in deriving energy estimates by which the nonlinear analysis is closed. 
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1.1. Equations and assumptions. We consider the general hyperbolic- 
parabolic system of conservation laws ([2]) in conserved variable C/, with 

-u G M, and v E M"~^, where, here and elsewhere, a denotes spectrum of 
a linearized operator or matrix. Here for simplicity, we have restricted to 
the case (as in standard gas dynamics and MHD) that the hyperbolic part 
(equation for u) consists of a single scalar equation. As in [MaZ3j . the results 
extend in straightforward fashion to the case u G M'^, /c > 1, with a{A^^) 
strictly positive or strictly negative. 

Following |MaZ4[ [Z3j . we assume that equations ([2]) can be written, al- 
ternatively, after a triangular change of coordinates 

w (u) 



(3) W := WiU) - , ^, ,,,, 



in the quasilinear, partially symmetric hyperbolic-parabolic form, 

(4) A'^Wt + AW, = (BW,), + G, 
where, defining W+ := W{U^), 

(Al) A{W+),A^,A^'^ are symmetric, A'^ block diagonal, A^ > Oq > 0, 
(A2) no eigenvector of A{A^)~^(W+) lies in the kernel of B{A^y^{W+), 

(A3) B=(^^ f\,b>9>0,a.ndG=(^}j with g{W,, W,) = 0{\W,\^). 

Along with the above structural assumptions, we make the following tech- 
nical hypotheses: 

{m)F,B,A^,A,B,W{-),g{;-)€C^. 

(HI) A^^ (scalar) is either strictly positive or strictly negative, that is, 
either A^^ > Oi > 0, oi A^^ < —6i < 0. (We shall call these cases the inflow 
case or the outflow case, correspondingly.) 

(H2) The eigenvalues of dF^^{U+) are real, distinct, and nonzero. 

(H3) Solution JJ is unique. 

Condition (HI) corresponds to noncharacteristicity, while (H2) is the condi- 
tion for the hyperbolicity of C/+. The assumptions (A1)-(A3) and (H0)-(H3) 
are satisfied for gas dynamics and MHD with van der Waals equation of 
state under inflow or outflow conditions; see discussions in |MaZ41 ICHNZi 
IGMWZ51IGMWZ6J . 

We also assume: 

(B) Dirichlet boundary conditions in M^-coordinates: 

(5) iw',w'')i0,t)=hit):={hi,h2m 



(7) 
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for the inflow case, and 

(6) w^^{0,t) = h{t) 

for the outflow case. 

This is sufficient for the main physical applications; the situation of more 
general, Neumann- and mixed-type boundary conditions on the parabolic 
variable v can be treated as discussed in [GMWZSl IGMWZ6] . 

Example 1.1. The main example we have in mind consists of laminar 
solutions {p,u,e){xi,t) of the compressible Navier-Stokes equations 

' dtp + div{pu) = 

dt{pu) + d\v{pu^u) + Vp = ep.Au + e{p + ?7)Vdivn 

dtipE) + div {{pE + p)u) = ekAT + epdiv{{u ■ V)n) 

-|- e{p + r])V{u ■ divn), 

X G M'^, on a half-space xi > 0, where p denotes density, u (^ W^ velocity, 

|2 

e specific internal energy, E = e + ^^- specific total energy, p = p{p, e) 
pressure, T = T{p,e) temperature, p > and \r]\ < p first and second 
coefficients of viscosity, k > the coefficient of heat conduction, and e > 
(typically small) the reciprocal of the Reynolds number, with no-slip suction- 
type boundary conditions on the velocity, 

Uj{0,X2,...,Xd) =0, j ^l and ui{0,X2, ■ ■ ■ ,Xd) = V{x) < 0, 

and prescribed temperature, T(0, X2, ■ ■ ■ , Xd) = Twaiiix). Under the standard 
assumptions pp, Te > 0, this can be seen to satisfy all of the hypotheses 
(A1)-(A3), (H0)-(H3); indeed these are satisfied also under much weaker 
van der Waals gas assumptions [ MaZ4[ [Z3l ICHNZl IGMWZ5[ [GMWZ6J . In 
particular, boundary-layer solutions are of noncharacteristic type, scaling as 
(/9, u, e) = (p, u, e){xi/e), with layer thickness ~ e as compared to the ~ ^/e 
thickness of the characteristic type found for an impermeable boundary. 

This corresponds to the situation of an airfoil with microscopic holes 
through which gas is pumped from the surrounding fiow, the microscopic 
suction imposing a fixed normal velocity while the macroscopic surface im- 
poses standard temperature conditions as in flow past a (nonporous) plate. 
This conflguration was suggested by Prandtl and tested experimentally by 
G.I. Taylor as a means to reduce drag by stabilizing laminar flow; see [S tlBraj . 
It was implemented in the NASA F-16XL experimental aircraft program in 
the 1990's with reported 25% reduction in drag at supersonic speeds [Brajo 
Possible mechanisms for this reduction are smaller thickness ~ e << ^/£ of 
noncharacteristic boundary layers as compared to characteristic type, and 
greater stability, delaying the transition from laminar to turbulent flow. In 



See also NASA site|http://www.dfrc.nasa.gov/Gallery/photo/F-16XL2/mdex.html| 
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particular, stability properties appear to be quite important for the under- 
standing of this phenomenon. For further discussion, including the related 
issues of matched asymptotic expansion, multi-dimensional effects, and more 
general boundary configurations, see |GMWZ5] . 

Example 1.2. For d?]), or the general ([2]), a large class of boundary-layer 
solutions, sufficient for the present purposes, may be generated as trunca- 
tions u^'^{x) := u{x — xq) of standing shock solutions 

(8) u = u{x), lim u{x) = u± 

on the whole line x G M, with boundary conditions /3/i(t) = u{0) (inflow) 
or Ph{t) = w {0) (outflow) chosen to match. However, there are also many 
other boundary-layer solutions not connected with any shock. For more gen- 
eral catalogs of boundary-layer solutions of ([7]), see, e.g., |MN1 ISZl ICHNZ^ 
l(4MWZ5j . 

Lemma 1.3 ( |MaZ3l [Z3l IGMWZSj l. Given (A1)-(A3) and (H0)-(H3), a 
standing wave solution ^ of ^, (B) satisfies 



(9) {d/dxfiU - U+) 

as X ^ -|-oo. 



< Ce"^^, A; = 0,...,4, 



Proof. As in the shock case |MaZ4t [Z3] . this follows by the observation that, 
under hypotheses (A1)-(A3) and (H0)-(H3), Uj^ is a hyperbolic rest point 
of the layer profile ODE. See also (GMWZl]. D 

1.2. Main results. Linearizing the equations ([2]), (B) about the boundary 
layer [/, we obtain the linearized equation 

(10) Ut = LU := -{AU)^ + (BU^)., 

where 

B:=B{U), AU := dF{U)U - {dB{U)U)U^, 
with boundary conditions (now expressed in [/-coordinates) 

(11) {dW/dif){Uo)U{0, t) = h{t) := (l'^ (t) 

for the inflow case, and 

(12) {dw^^/dU){Uo)U{0,t) = hit) 

for the outflow case, where {dW / dU){UQ) is constant and invertible, 

(13) {dw''/dU){Uo)=m{hi b2)m, 

(by (Al) and triangular structure ([3])) is constant with m G ]^("^i)><("^i) 
invertible, and h := h — h. 
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Definition 1.4. The boundary layer U is said to be linearly X ^ Y stable if, 
for some C > 0, the problem (jlOp with initial data Uq in X and homogeneous 
boundary data h = has a unique global solution U{-,t) such that |C/(-, i)|y < 
C|C/^q|x for all t; it is said to be linearly asymptotically X ^ Y stable if also 
\U{-,t)\Y ^0 as t^ oo. 

We define the following stability criterion, where D{X) described below, 
denotes the Evans function associated with the linearized operator L about 
the layer, an analytic function analogous to the characteristic polynomial of 
a finite-dimensional operator, whose zeroes away from the essential spectrum 
agree in location and multiplicity with the eigenvalues of L: 

(D) There exist no zeroes of D{-) in the nonstable half-plane ReX > 0. 

As discussed, e.g., in iR2l IMZTI iGMWZSi IGMWZ6J . under assumptions 
(H0)-(H3), this is equivalent to strong spectral stability, a{L) C {ReX < 
0}, (ii) transversality of C/ as a solution of the connection problem in the 
associated standing-wave ODE, and hyperbolic stability of an associated 
boundary value problem obtained by formal matched asymptotics. See 
|(;;MWZ5[ I(;MWZ6j for further discussions. 

Definition 1.5. The boundary layer U is said to be nonlinearly X ^ Y 
stable if, for each e > 0, the problem Q with initial data Uq sufficiently 
close to the profile U in \ ■ \x has a unique global solution U{-,t) such that 
\U{-,t) — U{-)\y < £ for all t; it is said to be nonlinearly asymptotically 
X ^ Y stable if also \U{-,t) — U{-)\y -^ as t ^ oo. We shall sometimes 
not explicitly define the norm X , speaking instead of stability or asymptotic 
stability in Y under perturbations satisfying specified smallness conditions. 

Our first main result is as follows. 

Theorem 1.6 (Linearized stability). Assume (A1)-(A3), (H0)-(H3), and 
(B) with \h{t)\ < Eo{l +t)~^. Let U be a boundary layer. Then linearized 
L^ n L^ ^>- L^ n U' stability, 1 < p < co, is equivalent to (D). In the case 
of stability, there holds also linearized asymptotic Li Ci L^ -^ L^ stability, 
p > 1, with rate 

(14) |?7(-,t)|z.P<C(l + t)~5(i-i/p)|t/o|iinLP + C7i?o(l + t)-^(i-i/p). 

To state the pointwise nonlinear stability result, we need some notations. 
Denoting by 

(15) tt]*" < a^ < • • • < a+ 

the eigenvalues of of the limiting convection matrix Aj^ := dF(U-^), define 

(16) e{x, t) := j; (1 + t)- V2e-|-4t|VMt^ 

a+>0 



(17) 
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Mx,t) := x{x,t) 5^ (1 + \x\+tr'/\l + \x- a+t\)-^/^, 



aj>0 



and 

(18) Mx, t) := (1 - x{x, t))(l + \x- a+t\ + t^/^'/^, 

where xix,t) = 1 for x G [0, a^t] and x{x,t) = otherwise and M > is a 
sufficiently large constant. 

For simplicity, we measure the boundary data by function 

2 

(19) Bn(.t):=Y,\{d/dtYh\' 

r=0 

for the outflow case, and 

4 2 

(20) Bh{t) ■.= ^\{d/dtYhi\^ + ^\{d/dtYh2\^ 

r=0 r=0 

for the inflow case. 

Then, our next result is as follows. 

Theorem 1.7 (Nonlinear stability). Assuming (A1)-(A3), (H0)-(H3), (B), 
and the linear stability condition (D), the profile U is nonlinearly asymp- 
totically stable in L^ n H'^, p > 1, with respect to perturbations Uq G H^, 
h ^ C^ in initial and boundary data satisfying 

||(l + |x|2)3/4[7o||^, <i^o and \BH{t)\ < E^{1 + \t\r^'^ 

for Eq sufficiently small. More precisely, 

\U(x,t)-U{x)\<CEo{9 + iPi + i^2){x,t), 
\U,ix,t)-U,{x)\ <CE^{e + i,i + ij2)[x,t), 

where U{x,t) denotes the solution of ([2]) with initial and boundary data 
U{x,0) = U{x) + Uo{x) and U{0,t) = Uq + h{t), yielding the sharp rates 

(22) \\U{x,t)-U{x)\\Lv<CEQ{l + t)-^'^^-^\ l<P<oo, 

(23) \\U{x,t)-U{x)\\LinH^<CE^{l + t)-\. 



Remark 1.8. By the one dimensional Sobolev embedding, from the hypoth- 
esis on Uq, we automatically assume that 

\\U4h^ < Eo, \Uo{x)\ + \U^{x)\ < Eo{l + |x|)-3/2. 

A crucial step in establishing Theorems 11.61 and 11.71 is to obtain pointwise 
bounds on the Green function G{x, t; y) of the linearized evolution equations 
pop (more properly speaking, a distribution), which we now describe. Let 
a^ ,j = 1, ....,n denote the eigenvalues of ^(+oo), and lt,r^ associated left 
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and right eigenvectors, respectively, normalized so that itr'^ = 6^. Eigenval- 
ues aj{x), and eigenvectors lj{x),rj{x) correspond to large-time convection 



rates and modes of propagation of the linearized model (jlOp . 
Define time-asymptotic, scalar diffusion rates 

(24) p+:={ljBrj) + , j = l,...,n, 
and local dissipation coefficient 

(25) r]^:=-D,{x) 
where 



D^{x) := Ai2b2^ 



A21 - A2262 ^h + &2 ^^1^* + b2dx{b2 ^h) (x) 



is an effective dissipation analogous to the effective diffusion predicted by 
formal, Chapman-Enskog expansion in the (dual) relaxation case, 

A^ ■.= An -Aub^^bi. 

Note that as a consequence of dissipativity, (A2), we obtain 

(26) r?+>0, /?+>0, for ah j. 

We also define modes of propagation for the reduced, hyperbolic part of 



system (flUj) as 

We define the Green function G{x, t; y) of the linearized evolution equa- 
tions (jlOp with homogeneous boundary conditions (more properly speaking, 
a distribution), by 

(i) (dt — Lx)G = in the distributional sense, for all x,y,t > 0; 

(ii) G{x, t; y) -^ 6{x — y) as t ^ 0; 

(iii) for all y, f > 0, I 7 * 7 j G(0, t;y) = (J where * = for the inflow 

case A^: > and * is arbitrary for the outflow case ^* < 0, noting that no 
boundary condition is needed to be prescribed on the hyperbolic part. 

By standard arguments as in |MaZ3j . we have the spectral resolution, or 
inverse Laplace transform formulae 

(28) e^V = —P.V. / e^*(A - Ly^fdX 
and 

-1 j-ri+ioo 

(29) G{x,t;y) = —P.V. e^'Gxix,y)dX 

^"^l Jrj—ioo 

for any large positive rj. 

We prove the following pointwise bounds on the Green function G{x, t; y). 
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Proposition 1.9. Under assumptions (A1)-(A3), (H0)-(H3), (B), and (D), 
we obtain 

(30) G{x,t-y)=H{x,t-y) + G{x,t-y), 

where 

and 

(32) 
\d2d^G{x,t;y)\ < Ce-''('^"2''+*) 



+ C'(t-(l"l + H)/2 ^ \^\^-v\y\ + |^|g-r,|x|)/ ^^-l/2g 



l/2^ix-y-a+tf/Mt 



k=l 
a+<0,a+>0 

< |a|, I7I < 1, /or some rj, C, M > 0, where indicator function X{|a+i|>|?/|> 

is 1 for |o^t| > \y\ and otherwise. 

Here, the averaged convection rate a*(x, t) in ([3T]) denotes the time-averages 
over [0, i] o/A^,(z) along backward characteristic paths z^, = z^,{x,t) defined 
by 

(33) —^ = A^{z^{x,t)), z^{t)=x. 

In all equations, a^, A^:, L^^^R^^ are as defined just above. 

1.3. Discussion and open problems. The stability of noncharacteristic 
boundary layers in gas dynamics has been treated using energy estimates 
in, e.g., [MNt IKNZJ IR3j . for both "compressive" boundary layers including 
the truncated shock-solutions dS]), and for "expansive" solutions analogous 
to rarefaction waves. However, in the case of compressive waves, these and 
most subsequent analyses were restricted to the small-amplitude case 

(34) 11^ — '"+IIli(ir+) sufficiently small. 

Examining this condition even for the special class ([8]) of truncated shock 
solutions, we find that it is extremely restrictive. 

For, consider the one-parameter family u^'''(x) = u{x — xq) of boundary- 
layers associated with a standing shock u of amplitude 5 := |ti+ — u_| << 1. 
By center manifold analysis [PeJ , n — u+ ~ 5e~ , hence 

11,1 „ II .-^-5^ \u+-u{0)\ 

^ ' \u+ — U-\ 

in fact measures relative amplitude with respect to the amplitude |u+ — U-\ 
of the background shock solution u. Thus, smallness condition (p4|) requires 
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that the boundary layer consist of a small, nearly-constant piece of the 
original shock. 

The present results, extending results of |YZ] in the strictly parabolic 
case, remove this restriction, allowing applications in principle to shocks 
of any amplitude. In particular, in combination with the spectral stability 
results obtained in [CHNZj by asymptotic Evans function analysis, they yield 
stability of noncharacteristic isentropic gas-dynamical layers of sufficiently 
large amplitude. Together with further, numerical, investigations of [CHNZj 
give strong evidence that in fact all noncharacteristic isentropic gas layers 
are spectrally stable, independent of amplitude, which would together with 
our results yield nonlinear stability. 

Spectral stability of full (nonisentropic) gas layers may be investigated 
numerically as for shocks in [HLyZl HLyZ2 , in both one- and multi- 



dimensions. However, analytical results of [SZ] show that in this case in- 
stability is possible, even for ideal gas equation of state. The numerical 
classification of stability for full gas dynamics, and the extension of our 
present nonlinear stability results to multi-dimensions, are two interesting 
direction for further investigation. 

Finally, we comment briefly on the difference between our analysis and 
the earlier analysis |YZ| carried out by similar techniques based on the Evans 
function and stationary phase estimates on the inverse Laplace transform 
formula. Our analysis is in the same spirit as, and borrows heavily from this 
earlier work. The main new issues are technical ones connected with the 
more singular high-frequency/short-time behavior of hyperbolic-parabolic 
equations as compared to the strictly parabolic equations considered in jYZj . 
In particular, linearized behavior in the u coordinate, U = {u,v), is essen- 
tially hyperbolic, governed for short times approximately by the principle 
part 

(35) vt + A,{x)v^ = 0, A, := {A^^^A^^ 

Thus, we may expect as in the whole-line analysis of hyperbolic-parabolic 
equations in |MaZ3j that the associated Green function contain a delta- 
function component transported along the hyperbolic characteristic 

dx/dt = ^*(x), 

with the difference that now we must consider also a possibly-complicated 
interaction with the boundary. 

A key point is that in fact this potential complication does not occur. 
For, in the special case occurring in continuum-mechanical systems [ Z3] 
that all hyperbolic signals either enter or leave the boundary, there is no 
such boundary interaction and no reflected signal. For example, in the 
simple scalar example (I35p . the Green function on the half- line with either 
homogeneous inflow {A^^ > 0) boundary condition f (0) = or outflow 
[A^"^ < 0) condition f (0) arbitrary, is by inspection exactly the whole-line 
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Green function 

9{x,t;y) = 6^-at{y)/A^{x) 

restricted to the half-line x,y > 0, where a is the average over [0, t] of 
A^:{z^{t)) along the backward characteristic path 

Indeed, comparing the description of the homogeneous boundary- value Green 
function in Proposition 11.91 with that of the whole-line Green function in 
|MaZ3] ■ we see that they are identical. However, to prove this simple obser- 
vation costs us considerable care in the high-frequency analysis. 

A further issue at the nonlinear level is to obtain nonlinear damping esti- 
mates using energy estimates as in [MaZ4] , which are somewhat complicated 
by the presence of a boundary. This is necessary to prevent a loss of deriva- 
tives in the nonlinear iteration. 

As in |YZ| . we get stability also with respect to perturbations in boundary 
data, something that was not accounted for in earlier works on long-time 
stability. We mention, finally, the works [GRl IMZTI ICJMWZSl IGMWZ6J in 
one- and multi-dimensions of a similar spirit but somewhat different techni- 
cal flavor on the related small viscosity problem- for example, e — > in ([7])- 
which establish that the Evans condition (or its multi-dimensional analog) 
is also sufficient for existence and stability of matched asymptotic solution 
as viscosity goes to zero. 

2. POINTWISE BOUNDS ON RESOLVENT KERNEL G\ 
In this section, we shall establish estimates on resolvent kernel Gx{x,y). 

2.1. Evans function framework. Before starting the analysis, we review 
the basic Evans function methods and gap/conjugation lemma. 

2.1.1. The gap/conjugation lemma. Consider a family of first order ODE 
systems on the half-line: 

W = A(x, X)W, Xen and a; > 0, 
^ '' M{X)W = 0, Xen and x = 0. 

These systems of ODEs should be considered as a generalized eigenvalue 
equation, with A representing frequency. We assume that the boundary 
matrix B is analytic in A and that the coefficient matrix A is analytic in A 
as a function from Q into L°°{x), C^ in x, and approaches exponentially to 
a limit A_|_(A) as x — > co, with uniform exponentially decay estimates 

(37) \{d/dxf{A - A+)| < Cie^^l^l/^^^ for x > 0, < k < K, 

Cj, 6 > 0, on compact subsets of il. Now we can state a refinement of 
the "Gap Lemma" of [GZl IKS| , relating solutions of the variable-coefficient 
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ODE to the solutions of its constant-coefficient limiting equations 

(38) Z' = K+{X)Z 

as X ^ +00. 

Lemma 2.1 (Conjugation Lemma |MZlj ). Under assumption (j37p . there 
exists locally to any given Aq € ^ a linear transformation P+(a;, A) = / + 
0+(a;, A) on X > Q, $+ analytic in A as functions from Jl to L°°[0,+oo), 
such that: 

(i) \P+\ and their inverses are uniformly hounded, with 

(39) \{d/d\y{d/dxfQ+\ < C{j)CiC2e~^^''^/^^ for x > 0, < k < K+1, 

j > 0, where < 9 < 1 is an arbitrary fixed parameter, and C > and the 
size of the neighborhood of definition depend only on 9, j, the modulus of the 
entries of A at Aq, and the modulus of continuity of A on some neighborhood 
o/Ao £n. 

(a) The change of coordinates W := P+Z reduces ()36p on x > to the 
asymptotic constant- coefficient equations ([381) . Equivalently, solutions of 
()36|) may be conveniently factorized as 



(40) w = ii + e+)z+, 

where Z+ are solutions of the constant- coefficient equations, and 0-|_ satisfy 
bounds. 

Proof. As described in [MaZ3j . for j = k = this is a straightforward 
corollary of the gap lemma as stated in [Z.3], applied to the "lifted" matrix- 
valued ODE 

P' = A+P - PA + (A - A+)P 

for the conjugating matrices P+. The x-derivative bounds < fc < i^ -|- 1 
then follow from the ODE and its first K derivatives. Finally, the A- 
derivative bounds follow from standard interior estimates for analytic func- 
tions. D 

Definition 2.2. Following [AG Jj . we define the domain of consistent split- 
ting for the ODE system W' = A{x, X)W as the (open) set of X such that the 
limiting matrix A_|_ is hyperbolic (has no center subspace) and the boundary 
matrix B is full rank, with dim 5+ = rankB. 

Lemma 2.3. On any simply connected subset of the domain of consistent 
splitting, there exist analytic bases {vi, . . . ,Vk}^ and {vk+i, ■ ■ ■ ,V]\f}~^ for 
the subspaces 5"+ and C/+ defined in Definition \2.Si 

Proof. By spectral separation of t/+, 5+, the associated (group) eigenprojec- 
tions are analytic. The existence of analytic bases then follows by a standard 
result of Kato; see pCat] . pp. 99-102. D 
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Corollary 2.4. By the Conjugation Lemma , on the dom,ain of consistent 
splitting, the stable manifold of solutions decaying as x ^ +00 of ([36]) is 

(41) S+:=span{P+vt,...,P+v+}, 

where W^. := P+v^ are analytic in A and C^^^ in x for A G C^ . 

2.1.2. Definition of the Evans Function. On any simply connected subset of 
the domain of consistent splitting, let W^ , . . . , W^ = P+v^ , . . . , P+V^ be 
the analytic basis described in Corollary 12.41 of the subspace S'^ of solutions 
W of (j36p satisfying the boundary condition 1^ ^ at +cxd. Then, the 
Evans function for the ODE systems W = A(x, X)W associated with this 
choice of limiting bases is defined as the k x k Gramian determinant 



L>(A) := det (j&W^, ..., MW+ 

= det (mP+v^, ..., MP.vt) 

V "^ / |a;=0,A 



(42) ^ ' 1"=°'' 



Remark 2.5. Note that D is independent of the choice of P^ as, by unique- 
ness of stable manifolds, the exterior products (minors) P+vf A • • • A P+v^ 
are uniquely determined by their behavior as x ^ +00. 

Proposition 2.6. Both the Evans function and the subspace S~^ are analytic 
on the entire simply connected subset of the domain of consistent splitting 
on which they are defined. Moreover, for A within this region, equation (|36]) 
admits a nontrivial solution W S L'^{x > 0) if and only if D{X) = 0. 

Proof. Analyticity follows by uniqueness, and local analyticity of P+, v'^ . 
Noting that the first P+v^ are a basis for the stable manifold of (j36p at 
X -^ +00, we find that the determinant of BP+f^ vanishes if and only if ]B(A) 
has nontrivial kernel on 5+ (A, 0), whence the second assertion follows. D 

Remark 2.7. In the case (as here) that the ODE system describes an eigen- 
value equation associated with an ordinary differential operator L, Proposi- 
tion {2^ implies that eigenvalues of L agree in location with zeroes of D. 
(Indeed, they agree also in multiplicity; see |GJH [GJ2] ; Lemma 6.1, [ZH]; 
or Proposition 6.15 of |MaZ3j . ) 

When ker B has an analytic basis v^,-^, . . . , t"^, for example, in the com- 
monly occurring case, as here, that B = constant, we have the following 
useful alternative formulation. This is the version that we will use in our 
analysis of the Green function and Resolvent kernel. 

Proposition 2.8. Let v^,-^, . . . ,v^j^ be an analytic basis o/kerB, normal- 
ized so that det (B*, w^^^, . . . f^) = 1. Then, the solutions W^ of (f36]l 
determined by initial data W^{X,0) = v^ are analytic in A and C^^^ in x, 
and 

(43) Z)(A):= det (w^+,...,l^+,l^fcV>- ••><), , ,■ 
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Proof. Analyticity/smoothness follow by analytic/smooth dependence on 
initial data/parameters. By the chosen normalization, and standard prop- 
erties of Grammian determinants, 

D{X) = det(w^,...,W^,vl_,„...,v%) 

\ / |a;=0,A 

yielding dM]). D 

2.1.3. The tracking/reduction lemma. Next, consider a family of systems 
W' = A{x,p,e)W, pGP, eGM+ and x > 0, 
B{p,e)W = 0, Xen and x = 

parametrized by p, e, with e — > 0. The main example we have in mind is 
()36p with p = A/|A| and e := lAj^"*^, in the high-frequency regime |A| -^ oo. 
We assume further that by some coordinate change we can arrange that 

(45, A^(1^.?>e. 

with 

(46) |e| < 6{e), K(M+ - M_) > 2r]{e) + a^(x), 
I|q^IIli(r+) uniformly bounded for all e sufficiently small, and 

(47) {6/r]){e)^0 as e ^ 0, 

where ^{Q) := (1/2) (Q + Q*) denotes the symmetric part of a matrix Q. 

Then, we have the following analog of Lemma 12.1.11 asserting that the 
approximately block-diagonalized equations ()44p may be converted by a 
smooth coordinate transformation 

92 / 



,2 T I ^ J as e ^ 



to exactly diagonalized form with the same leading part M. 

Lemma 2.9 ( [MaZSj ). Consider a system (j45p . with F = and 5/r] -^ 
as e -^ 0. Then, (i) for all < e < eo, there exist (unique) linear 
transformations ^l{z,p) and <l>|(z,p), possessing the same regularity with 
respect to the various parameters z, p, e as do coefficients M± and 0, for 
which the graphs {(Zi, <I>|Zi)} and {{^\Z2, Z2)} are invariant under the 
flow of (I45p . and satisfying 

l^il, 1^2! < C6{e)/ri{e) for all z. 

In particular, (ii) the suhspace -E_ of data at z = for which the solution 
decays as z ^ +00, given by span {{^1(0, p)v,v)}, converges as e ^ to 
-E_ := span{(0, u)}. 

Proof. Standard contraction mapping argument carried out on the "lifted" 
equations governing the flow of the conjugating matrices <l>j; see Appendix 
C, |MaZ3j . 
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Remark 2.10. In practice, we usually have a^ = 0, as can be obtained 
in general by a change of coordinates multiplying the first coordinate by 
exponential weight eJ '^^'^^. 

2.2. Construction of the resolvent kernel. In this section we construct 
the explicit form of the resolvent kernel, which is nothing more than the 
Green function Gx{x, y) associated with the elliptic operator (L — XI), where 



(48) (L - XI)Gx{;y) = 6yl, ^-^* Q Gx{0,y) ^ ^^ 

where * = for the inflow case and is arbitrary for the outflow case. 

Let A be the region of consistent splitting for L. It is a standard fact (see, 
e.g., [He]) that the resolvent {L — XI)~^ and the Green function G\{x,y) are 
meromorphic in A on A, with isolated poles of finite order. 

Writing the associated eigenvalue equation LU — XU = in the form 
of a first-order system ()36p as follows: W := {u,v,z) G C^"~^ with z := 
biu' + b2v' , and 

u' = A-\-Aub^^z - {A[, + X)u - A[^v), 

(49) v' = h^^z - b^^biu', 



z 



{A21 - A2262 ^bi)u' + A2262 ^z + A'^iu + {A'22 + X)v. 



2.2.1. Domain of consistent splitting. Define 

(50) A:=nA+, j = l,2,...,n 

where A^ denote the open sets bounded on the left by the algebraic curves 
A^(^) determined by the eigenvalues of the symbols —^^Bj^ — i^Aj^ of the 
limiting constant-coefficient operators 

(51) L+w := B+w" - A+w' 

as ^ is varied along the real axis. The curves A^ comprise the essential 
spectrum of operators -L+. 

Lemma 2.11 ( |MaZ3j ). The set A is equal to the component containing real 
+CXD of the domain of consistent splitting for (j49p . Moreover, if assumption 
(H3) holds, then 

(52) Ac{A : 3[?eA> -r/|9mA|/(l + |9mA|), r? > 0. 

2.2.2. Basic construction. We first recall the following duality relation de- 
rived for the degenerate viscosity case in |MaZ3j . 

Lemma 2.12 f [ZHl IMaZ3J 'l . The function W = {U, Z) is a solution of (09]) 
if and only if W*SW = constant for any solution W = (U, Z) of the adjoint 
eigenvalue equation, where 

-An -A12 0' 

(53) S = I -A21 -A22 Ir 

-b^^bi -Ir 
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and 

(54) Z = {hiM)U\ Z = {0,bl)U'. 
For future reference, we note the representation 

/ -A-^ A-^Au 

(55) 5-1 = b^%A-^ -62-^61^-1^12 - Ir 

V -AA-^ Ir -A22 + AA-^Au 

where A := A21 — A226-161, A^, := An — ^1262^ 61, obtained by direct 
computation in [MaZ3] . 
Denote by 

(56) c^° = (<A^+i(x;A),---,<A°+,(x;A)), 

(57) <^+ = i<Piix;X),... ,cl>tix;X) = iP+vtr-- ,P+vt), 
and 

(58) $ = ($+,$°), 

the matrices whose columns span the subspaces of solutions of ([36]) that, 
respectively, decay at a; = +00, and satisfy the prescribed boundary con- 
ditions at a; = 0, denoting (analytically chosen) complementary subspaces 
by 

(59) M/0 = (V'?(x;A),---,V'0(x;A)), 

(60) *+ = (V'^_,,(x;A),...,V++,(x;A)) 
and 

(61) ^ = (^°,^+). 

As described in the previous subsection, eigenfunctions decaying at +00 
and satisfying the prescribed boundary conditions at occur precisely when 
the subspaces span$^ and span<l>"^ intersect, i.e., at zeros of the Evans 
function defined in (14311: 



(62) Di(A):=det(cI.o,<I>+)|,=o- 

Define the solution operator from y to x of (L — A)[/ = 0, denoted by 
J^S'-^^, as 

Ty^- = ^x,X)<^-\y,X) 

and the projections nlj,!!^ on the stable manifolds at 0,-|-oo as 

n+ = ($+(2/) 0)'^-\y), nO = (0 $0(y))cl>-i(y). 

With these preparations, the construction of the Resolvent kernel goes 
exactly as in the construction performed in |ZHl IMaZ3] on the whole line. 
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Lemma 2.13. We have the the representation 

,.„^ r( ^iiIn,0)Ty-^U+S-Hy){In,Or\ for x > y, 

^ ' '^^'^^ \-(/„,0)^^--^05-l(y)(/n,0)*^ for x<y. 

Moreover, on any compact subset K of p{L) f] A, 

(64) \Gx{x,y)\<Ce^^^-y\, 
where C > and r/ > depend only on K,L. 

We define also the dual subspaces of solutions of (L* — X*)W = 0. We 
denote growing solutions 

(65) ^o = {4>0(x;X) ■■■ 4>lix;X)), 

(66) l>+ = (<A++i(x;A) ••• <A++.(^;A)), 
$ := ("J)*^, $+) and decaying solutions 

(67) §° = (v;?(x;A) ••• ^+(x;A)), 

(68) *+ = (^+^,(x;A) ••• i^t+rix;X),) 
and ^ := (^'^,^'''), satisfying the relations 

{M)U{^^)o,+ ^ I. 

Then, we have 
Proposition 2.14. The resolvent kernel may alternatively be expressed as 

,.q. ^f f(I„,0)cl>+(x;A)M+(A)§o*(y;A)(^,0)*'^ x > y, 

^ ^ '^ '^' \-(I.,0)$0(x;A)MO(A)§+*(y;A)(/„,0)*- x<y, 
where 

(70) M(A) :=diag(M+(A),MO(A)) = '^-\z;X)S-^{z)^-^*{z;X). 

From Proposition 12.141 we obtain the following scattering decomposition, 
generalizing the Fourier transform representation in the constant-coefficient 
case 

Corollary 2.15. On An p{L), 

(71) Gx{x, y) = Y, d%'Pt{x; A)V5fc+(y; A)* + J^ c/>+(x; A)(^+(y; A)* 

j,k k 

for < y < X, and 

(72) Gx{x,y) = Y,d%{X)(Pt{x;X)i^+{y;Xr + J24^+{x;X)i^-l{y;Xr 

j,k k 

for < X < y, where dA (A) = 0{X~ ) are scalar meromorphic functions 
with pole of order K less than or equal to the order to which the Evans 
function D{X) vanishes at X = (note that K = under assumption (D)). 
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Proof. Matrix manipulation of expression ()70p . Kramer's rule, and the def- 
inition of the Evans function; see |MaZ3| . D 

Remark 2.16. In the constant- coefficient case, with a choice of common 
bases $^'+ = <I)+'0 at 0, -|-cxd, the above representation ()2.15p reduces to the 
simple formula 

2.3. High frequency estimates. We now turn to the crucial estimation 
of the resolvent kernel in the high-frequency regime |A| -^ +oo, following 
the general approach of |MaZ3j . Define sectors 

(74) Qp := {A : ^eX > -0i|9mA| + 62}, 9j > 0. 
and 

(75) := {A : -??i < 5ReA} 

with r]i sufficiently small such that O,\B(0,r) is compactly contained in the 
set of consistent splitting A, for some small r to be chosen later. Then, we 
have the following crucial result analogous to the estimates on the whole 
line performed in [MaZ3] . 

Proposition 2.17. Assume that (H0)-(H3) hold. Then for any r > and 

"Hi — ^il'') > chosen sufficiently small such that fi \ B{0,r) C A n p{L). 
Moreover for R > sufficiently large, the following decomposition holds on 
n\B{0,R): 

(76) Gx{x,y) = H^{x,y) + Px{x,y) + ef (x,y) + ef (x,y), 
where 

and 
(78) 
where 

(79) «.(...) =c.(.,.)=^'->»'^;'::;:::;:;::'-<^-='' ^>''- 

[X{A,<o}e ■'y '' ' b^(x,y) x<y, 
with 

(80) b, := j:(-v*/^*)i-)d- = 0(e-^l"-^l), 
due to (j26p . and 

(81) D^{x,y;X) = o(e-e(i+KeA)l.-,l ^ g-eiA|V2|,_,|^^ 



ef (x, y) = \-^Bxix, y; A) + X-\x - y)Cx{x, y; A), 
@^{x,y) = X-^Dx{x,y;X) 
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for some uniform 9 > independent ofx, y, z, each described term separately 
analytic in A, and P\ is analytic in X on a (larger) sector Qp as in (I74p . with 
6i sufficiently small, and 62 sufficiently large, satisfying uniform bounds 

(82) {d/dx)''{d/dyfPx{x,y) = 0(|A|(l°l+l^l-i)/2)g-e|A|i/2|^_^|^ ^ ^ ^^ 

for \a\ + |/3| < 2 and < |a|, |/3| < 1. 

Likewise, the following derivative bounds also hold: 

{d/dx)@x{x, y) ={bI{x, y; X) + {x - y)Cl{x, y; A)) + X'^ [bI{x, y; X) 
+ {x-y)Cl{x,y;X) + {x - yfDl{x,y; X)^ 

+ X-^/^E^{x,y;X) 
and 

{d/dy)@x{x, y) = (i?°(x, y; A) + (x - y)Cl{x, y; A)) + A-^ (i?i(x, y; A) 

+ (x - y)C^(x, y; A) + (x - yfDl{x, y; A)^ 
+ X-''/^Ey{x,y;X) 

where Bn, Cp, and Di satisfy bounds of the form ([79]) . and Ep satisfies a 
bound of the form (j8ip . 

Proof. We shall follow closely the argument in |MaZ3j . with the new feature 
of boundary treatments, or estimates of $*^,^''. Writing the associated 
eigenvalue equation LU — XU = in the form of a first-order system as 
follows: W := {u,v,z) G C^"^-*^ with z := biu' + b2v' , and 

u' = A-\-Ai2b2^z - (A'li + X)u - A[^v), 

(83) v' = b^^z - b:^^biu', 

z' = {A21 - A22b2^bi)u' + A22b2^z + A'^iu + {A'22 + X)v 
or 

(84) W' = AW. 

Recall from Lemma 12.131 that we have the the representation 

.o,^ ^ . .^j{In,0)J^'w''U^{y)S-\y){In,Or, for x > y, 

^ ' '^ '^^ \-(/„,0)^^-^^^o^(y)5-l(y)(/„,0)*^ for x<y. 

We shall find it more convenient to use the "local" coordinates u := 
A^:U,v := biu + b2V. yielding from (P9]) : 

Ux = -XA^^U - {Ai2b2^v)::i: 

(86) (v,), = [{{A21 - A22b2\ + b2d,ib2^bi))A-'u), 

+ ((^22 + dx{b2)b^^)i)x + Xb^^biA-^u + Xb^^i 
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Following standard procedure (e.g., |AGJ1 iGZt IZHl [MaZ3] ). performing 
the rescaling 

(87) x:=\X\x, A:=A/|A|, 

and changing coordinates W >-^ Y = QW, where 

Y = {u, V, VxY'' = {A^u, biu + b2V, {bm + 62w)x•)*^ 

/A, 

Q= \ bi 62 

\\xr'dxbi \\\~'dxb2 lAi-i/,,, 



and 










/ 


A-' 





(90) 


Q- = 


-b2%A-' 


b^' 






-\X\b2Mb^%)A-' -\xmb2)b2' |A|/r, 

we obtain the first order equations 

(9f) Y' = Aix,\X\-^)Y, Y:=iu,v,vT, '-=9, 

where 

(92) A{x, |A|-i) = Ao{x) + \X\-'Aiix) + 0{\X\-^), 

with 



(93) 



An(x) 



Aiix) 



XA-^ 




-Ai2b^^\ 

'0 ) 







Xd^A-^ 


-d.iAub^^) 



A6^^ 







(94) 



d^ := A21 - A2262 ^61 - 62 ^biA^ + b2dx{b2 %), 
e* := A22 + d^A-^Au + 5^(62)- 



We will carry out the details of the lower-order estimates in Proposition 
12. 17^ leaving high-order estimates and derivative bounds as brief remarks at 
the end. First, observe that the representation (f85|l becomes 
(95) 

r( N^ /(^n'0)2"'•^^"^J(y)2'^-l(y)(/„,0)*^ for x>y, 

'^"^'^^ \-(/„,0)Q-l.Fr^■^o.(y)Q^-^y)(^,0)*^ for x<y 

where Ily and T^^ denote projections and flows in y— coordinates. 
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2.3.1. Initial diagonalization. Applying the formal iterative diagonalization 
procedure described in |MaZ31 Proposition 3.12], one obtains the approxi- 
mately block-diagonalized system 

(96) Z' = D{x, \X\~^)Z, TZ := Y, D := T-^AT, 

(97) Tix, \Xr') = Toii) + \X\-^T,ii) + ■■■ + \X\-'Ts{x) 

(98) D{x, \X\-') = Do{x) + \X\-'Di{i) + ■■■ + D^{x)\X\-^ + 0{\X\-^), 

where without loss of generality (since Tq is uniquely determined up to a 
constant linear coordinate change) 

/l -X-^A^Ai2b2^\ (l X~^A^Ai2h2^' 

(99) To := /, 1 Tq-i = I /, 

\0 /,. / 

and 



Ir \0 Ir 




(100) Dq:=\ / J , £> 



'1 ■-- 




with r]^ as defined in ()25p : see Proposition 3.12 |MaZ3j . (Here, the simple 
block upper-triangular form of ^o has been used to deduce the above simple 
form of Do, Di.) 

2.3.2. The parabolic block. At this point, we have approximately diagonal- 
ized our system into a 1 x 1 hyperbolic block with eigenvalue fl = —X/A^, of 
^0; and a 2r X 2r parabolic block 

(101) 4 = NZp 
with 

(102) A':=(S t)+l^|-'G"-. ")+°(W-')- 

Balancing this matrix A^ by transformations B := diag{Ir, \X\^'^''^Ir} we 
get 

(103) M := B-^NB = lAr^/^Mi + 0{\X\-^), Mi := L^_i M 



Observe that cf{Mi) = ity cr(A62 ) has a uniform spectral gap of order 
one. Thus, there is a well-conditioned transformation S = S{Mi) depending 
continuously on Mi such that 

(104) Ml := S-^MiS = diag{ Af",M+} 

with M-^ uniformly positive/negative definite, respectively. Applying this 
coordinate change, and noting that the "dynamic error" S^^dxS is of order 
dxMi = 0(1 A|~"^), we obtain the formal expansion 

(105) M{x,\X\-') = \X\-'/^dieig{M{,M^} + (D{\X\-'). 
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Finally, on sector i7p, blocks |A|~^'^Af-|^ are exponentially separated to 
order |A|~^'^. Thus, by the reduction lemma, Lemma |2.9[ there is a further 
transformation S := l2r + 0(1 A|^^'^) converting M to the fully diagonalized 
form 

M(x,|A|-i):=|Ari/25-i(Mi + 0(|A|-i/2)^5 

= 0(|Ari/2)diag{Mf,M+} 

where Mj^ = Mj^ + 0(1 A|~^'^) are still uniformly positive/negative definite. 
In summary, changing coordinates 

(106) BSSZp = Zp, 
(dnH) yields 

(107) ^; = 0(|ArV2) (^Mf ^^^^Zp + 0{\Xr'/') 

Therefore the transformation 

(108) T:=(To + |A|-iri)(j ^°^ 

converts equations (|91|) to the following: 

c = -(aa;i + \x\-%a:')c + odAr^) 
^ ^ ,'i = iArV2Mf^^ + o(iAr3/^) 

by relation 

(110) TZ = Y, Z = {C,P-,P+T- 

Then, we have the the representation 
(111) 

, f(I„,0)Q-lr^r^^+(y)r-lQS-l(y)(/„,0)*^ for x>y, 

'^"^'^^ \-(I„,0)Q~lr^r'^|(y)T-lQ5-l(y)(/n,0)*^ for x<y, 

thanks to the fact that 

(112) ^^""" = T^|"""T-\ n+ = rn+T-i. 

Computing, we have 

/I |Ari/2 |A|-i/2\ A A-M,Ai262-^^ 

T= 0(1) 0(1) T-i= 0(1) |A|i/2 

\0 |Ari/2 |A|-V2/ \o 0(1) |A|i/2 

and 

(1^^) (^-0)2" = U4U- 6,- o) 

ai4) ri oiP-iT-T ^*' o(IAr^/') o(|Ari/^)A 
(114) (/„,0)Q r-^_^_,^^^_, ^^^^ ^^^^ j 
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and 

(115) QS~\ln,0Y' 



/_ l + |A|-i 0(|A|-i) 

(116) T-'QS-\ln,Of=\o{\\\-'/^) 0(|ArV2) 

VO(|A|-V2) o(|A|-V2) 

Therefore now we are ready to estimate ^z~^^^% ^^^ -^^^^z- 

2.3.3. Estimates on projections and solution operators. We shall give esti- 
mates on the projections: 

(117) n+ = ($+, o)($+, $°)-i, U% = (0, $«)($+, $0)-i 

and the solution operators: 

(118) .^r^ = (cl>+(x), cl>0(x))(cl>+(y), $0(y))-i. 
First, let (I)P+/\I/P+ be the decaying/growing basis solutions of 

(119) p'_ = |Ar^/2^f p_ and p'^ = \X\-^/^M+p+ 
and (j) ~^ ftp ~^ be the decaying/growing basis solutions of 

(120) C = -(A^-i + \X\'%A:')C. 

Lemma 2.18. [Inflow case] For the inflow case A* > 0, we obtain 

(I -|A|-V20?'+e(A)^P+'^^ 

(121) n+ = ir -$p+£;(A)^p^ 

Vo 

/O |A|-V20^+e(A)^P+-i^ 

(122) n| = $P+S(A)*P+"^ 

Vo Ir 

wit/i hounded functions e{X),E{X), and 

(<t)^+{x)^^+{y)-^ 

(123) .F|-"^ = $P+(x)$P+(y)-i 

\ ^P+(x)^P+(y)~\ 

Proof. We have the decaying basis solution in 2^-coordinates of the first 
order equations (|109p 

h'^+ 

(124) $+= ( $P+ I +0(|Ar^) 
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Since $"*" and ^^ (exactly ^P"*") form a basis solution, we can write 

(<i^^+\ / \ / 

(125) $°(x) = e(A) + <^P+{x)E{\) + 

\ / V / \^P+{x)F{\)^ 

Now since {ij/- }j forms a basis, we can take {ij/- (0)} to be the analytic 
basis for y at x = 0. Also as we recall that Z = T^^Y, we compute 

/ \ / 0(1) 

(126) ^1-0 = '^-' { = |A|^/Vr(0) 

This and (fT25]l yield 

/ e(A)0'^+(x) \ 

(127) -^"(x) = \X\^/^^P+ix)E{X) +0{\X\-^/^) 

\ |A|i/2$P+(x) / 

where 

(128) E{X) = {Ei{X),...,Er{X)f\ i?,(A) = ^f (0,A)$P+(0,A)-i 

and e(A),£'j(A) G M^ are bounded functions in A. Therefore computing, we 
get 

U^+'^ -|A|-V2e(A)^p+-i' 

(129) ($+,$0)-i= $P+"^ -E{X)^P+-^ 

\ |A|-V2^p+-i 

and hence straightforward computations give the lemma. D 

Lemma 2.19. [Outflow case] For the outflow case A* < 0, we obtain 
/O \ /^ ^ 

(130) n+ = /,. -$p+^^p+-^ , n| = ^p+ew+~^ | , 

\0 / \0 /^ 

where E{X) is a bounded function in X determined below. Moreover, 

/ij^+{x)ij^+{y)~^ 

(131) J^^'-" = $P+(x)$P+(y)-i 

\ W+{x)^P+{y)-^^ 

Proof. Similarly, we have $"*" = ^^"^ and ^^ = {(1)^^ , $pO) where we can write 

(132) $"(x) = I ^P+{x)E{X) I +e(A) ( | + 

\^P+{x)F{X)^ 
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As before, using the form of the hnearized boundary conditions (jT2|l . we 
can take 

/ \ / 0(1) 

(133) <^^o =T-M = |A|VVf (0) 

and thus 

(134) $P°(x) = |A|V2$p+(x)^(A) 

V |A|V2$P+(x) / 

with bounded functions e(A) and ^^(A) = V^~^(0, A)$p+(0, A)-^ 
Similarly, we take 

and thus 

(135) (/>^°(a;) 

Putting together and computing, we obtain 

/ i^+ e(A)^'^+ \ 

(136) ($+,$°)=Up+ |A|V2$p+s(A) 

V |A|V2^p+ y 

and 

/ $P+"^ -^(A)^P+~^ 

(137) ($+,$0)-i=L^+-i -|A|-V2e(A)^P+ 

\ |A|-V2^p+-i 

Direct computations yield the lemma. D 

2.3.4. Estimates on G\: Inflow case A^ > 0. Now we are ready to combine 
all above estimates to give the bounds on resolvent kernel Gx. We shall 
work in detail for the case x > y. Similar estimates can be easily obtained 
for X < y. First decompose the projection as Ili = 11^^ + 11?, where 




(138) 



TTfe+ 



TTP+ 




z 

IN 
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Hence 
H^{x,y) = (/„,0)Q-iT^r^ni+(y)T-iQ5-i(y)(/„,0) 



tr 






^b^'hA-'{x) 0^ 
= cp''+{x)^^+iy)-'RM: + 0{\Xr')cP^+{x)cP^+iyr\ 

recalling that 0^+(x)0^^(y)^^ is the solution operator of hyperbolic equa- 
tion in (I120p and thus satisfies 
(139) 

(/)^+(2;)(/)^+(y)-i = e4"(-iM*-|A|-i»7.M*)(^)d2 ^ j;{-x/A,-ri,/A,){z)dz _ 

At the same time, computing Px{x,y), we obtain 

PA(:r,y) = {In,0)Q-'TJ^r"^''+{y)T-'QS-\y){In,Oy' 

= 0{\X\~^/^)^P+{x)^P+{y)-^ 

recalling that ^'P~^ (x)^^^ {y)^^ is the (stable) solution operator of parabolic 
equation (J119p . with M^ uniformly negative definite, and thus we have an 
obvious estimate 

(140) |$P+(x)$P+(y)-^| < C7e-^l^l"'^'(^-*) < Ce^^l^l'^'^^-^). 

We therefore obtain 

(141) P,(x,y)=0(|A|-V2)g^e|A|V2(,_,)^ 

2.3.5. Estimates on G\: Outflow case A^ < 0. Again as above, we shall 
work in detail for the case x > y. Similar estimates can be easily obtained 
for X < y. Estimates in Lemma l2. 191 yield 

/O 

(142) J^|~"^n+(y) = <^P+{x)<^P+{y)-^ -$P+(x)E(A)^P+(y)-i 

\0 

where $P+(a;)^(A)^P+(y)-i < C^P+{x)^P-^{y)-K Observe that n|+ = 0. 
Therefore, H\{x,y) = and 

Py^{x,y) = (/„,0)Q-iT.Fr"n|+(y)T-iQ5-i(y)(I„,0)*'- 



< c|ArV2e-^l^l'^'(-?') 

We thus complete the proof of estimates Hx and of Px appearing in Propo- 
sition EUl 
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2.3.6. Derivative estimates. Derivative estimates now follow in a straight- 
forward fashion, by differentiation of (jllip . noting from the approximately 
decoupled equations that differentiation of the flow brings down a factor (to 
absorbable error) of A in hyperbolic modes, A^' ^ in parabolic modes. This 
completes the proof of Proposition 12. 171 D 

2.4. Low frequency estimates. Our goal in this section is the estimation 
of the resolvent kernel in the critical regime |A| — > 0, i.e., the large time 
behavior of the Green function G, or global behavior in space and time. We 
are basically following the same treatment as that carried out for viscous 
shock waves of strictly parabolic conservation laws in [Zm IMaZ3j ; we refer 
to those references for details. In the low frequency case the behavior is 
essentially governed by the limiting far-field equation 

(143) Ut = L+U := -A+U^ + B+U,, 

Lemma 2.20 ([MaZ3]). Assuming (H0)-(H3), for |A| sufficiently small, the 
eigenvalue equation {Lj^ — A)VF = associated with the limiting, constant- 
coefficient operator L^, considered as a first- order system W = A^W , W = 
{u,v,v'), has a basis of 2n — 1 solutions W^ = e^+^^'^Vj{X), consisting of 
n — 1 "fast" modes (not necessarily eigenmodes) 

(144) |e*+(^)^yj| < Ce^^l^'l, 9 > 0, 
and n analytic "slow" (eigen-)modes 

(145) /^ti+,(A) := -A/a+ + X^dpaf + 0{X% 
V;ti+,(A):=r; + 0(A), 

where a'j' , l^ , r^,f3^ are defined as in Proposition \1.9l The same is true 
for the adjoint eigenvalue equation 

(L+ - xyz = 0, 

i.e, it has a basis of solutions W^ = e~ +^ '^Vj{X) with n — 1 analytic "fast" 
modes 

(146) |e-^+(^)^Fj| <C7e-^l^l, 6 > 0, 
and n analytic "slow" (eigen-)modes 

(147) V;+_i+/A) = /+ + 0(A). 

Proof. Standard matrix perturbation theory; see [MaZ3] . Appendix B. D 
Also we recall from the representation of G\ in Corollary 12.151 
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Proposition 2.21. Assuming (H0)-(H3), let K be the order of the pole of 
G\ at X = and r be sufficiently small that there are no other poles in 
B(0,r). Then for A € f^g such that |A| < r and we have 

(148) Gx{x,y) = J^dJ (A)<A+(x)^+(y) + J],/.+ (x)0+(y), 

j,k k 

for X > y > 0, and 

(149) Gx{x,y) = Y,d%{X)(P+{x)i^+{y) + J2i^^i^)i'Uyl 

j,k k 

for < X < y, where dA (A) = 0{\^ ) are scalar meromorphic functions, 
moreover K < order of vanishing of the Evans function D{X) at X = 0. 

Proof. See |ZH1 Proposition 7.1] for tiie first statement and theorem 6.3 for 
the second statement hnking order K of the pole to multiphcity of the zero 
of the Evans Function. D 

Our main result of this section is then: 

Proposition 2.22. Assume (H0)-(H3) and (D). Then, forr > sufficiently 
small, the resolvent kernel G\ associated with the linearized evolution equa- 
tion (|143p satisfies, for < y < x: 
(150) 

\d2d^Gx{x,t;y)\ 

<C(|Ar + e-^l^l)(|Ar + e-^l^l)( ^ |e(-A/a++AX+/4')(-?/)| 



+ E 



a+>0 

I (-A/a++A2/3+/a+')x+{A/a+-A2/3+/af)j/ 



a+<0,a+>0 

< |a|, I7I < 1, > 0, with similar bounds for < x < y. Moreover, each 
term in the summation on the righthand side of ()150p bounds a separately 
analytic function. 

Proof. By condition (D), D(X) does not vanish on Re{X) > 0, hence, by 
continuity, on |A| < r. Thus, according to Proposition 12. 211 all \djk{X)\ are 
uniformly bounded on |A| < r, and thus it is enough to find estimates for 
fast and slow modes i;^^, 0^, tp^ and t/j^ . By applying Lemma 12.201 and 
using (f57|) we find: 

and similarly for 0^, i/j^ and ^p^ . Now using ([39]) and the fact, by Lemma 

[2:201 that e^j(^)^ is of order e~'^^^''t+^^'^t/"'t +C(A'))^ for slow modes and 
order e~ '^' for fast modes, so by substituting this and corresponding dual 
estimates in (I15ip and grouping terms, we obtain the result. 

D 
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3. POINTWISE BOUNDS ON GrEEN FUNCTION G{x,t;y) 

In this section, we prove the pointwise bounds on the Green function G 
fohowing the general approach of [MaZ3] in the whole-hne, shock, case. Our 
starting point is the representation 

1 r-q+ico 

(152) G{x,t;y) = —P.V. e^'Gx{x,y)dX 

^T^i Jrj—iQO 

where r] is any sufficiently large positive real number. 

Case I. \x — y\/t large. We first treat the simple case that \x — y\/t > S, 
S sufficiently large. Fixing x, y, t, set A = r/ + i^, for 77 > sufficiently large. 
Applying Proposition 12.171 we obtain the decomposition 

G(x, t; y) = —P.V / e^* //^ + ef + Pa + ef (x, y)d\ 

= : I + 11 + III + IV. 

For definiteness considering the inflow case A* > and taking x > y, we 
estimate each term in turn. 

Term I. Computing, 

(•ri+ioo 



1 nri^tOQ 

I=—P.V e^'Hx{x,y)dX 

^T^l J rj—ioo 



2vr * 7_oo 



2TT Jy 

zvr 
where a^, is defined as in Proposition II. 9[ Noting that a* > inf^. ^*(x) > 
and r/+ > 0, we get g" ■/'«('?* /^*)(^)'^^ = ©(e^^^^-J^)) and thus 

(153) / = 0{e-'')6^.-a,t{y), 

vanishing for |x — y\/t large. 

Term II. Similar calculations show that the "hyperbolic error term" // also 
vanishes. For example, the term e A~^i?(x,y;A) contributes 

27r 7_oo 

The integral though not absolutely convergent, is integrable and uniformly 
bounded as a principal value integral, for all real rj bounded away from zero. 
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by explicit computation. On the other hand 

^V(t-I;i/A,(z)dz) ^ ^,j{t~\x-y\/min,A4z)) < ^^1(1- S / min, A^z)) _^ q^ 

as r/ — > +00, for S sufficiently large. Thus, we find that the above integral 
term goes to zero. Likewise, the result applies for the term of e^*C(a;,y; A), 

since {x — y)e~-'y ^'^*' *'^^' ^ < C{x — y)e~^^^~y' is also bounded. Thus, each 
term of II vanishes as r/ — > +00. 

Term III . The parabolic term /// may be treated exactly as in the strictly 
parabolic case |ZHj . Precisely, we may first deform the contour in the prin- 
ciple value integral to 



Jv^i 



(154) / e'^'Px{x,y)d\ 

lUra 

where Fi := dB{0, R) D Clp and r2 := di^p \ B(0, R), recalling the parabolic 
sector Q,p defined in (i71|) . Choose 

(155) «:=t_^, R:=9a^ 



where 9 is as in (|82|) . Note that the intersection of T with the real axis is 
Amin = R = Qo? . By the large |A| estimates of Proposition [2T71 we have for 
all A G Fi U Ps that 

|P,(x,y)|<C|A|-VV^I^I^^^I-< 
Further, we have 

(156 - V / ;, 

Re\ < ReXo - 7?(|/mA| - |/mAo|), A G F2 

for R sufficiently large, where oj is the argument of A and Aq and Ag are the 
two points of intersection of Fi and F2, for some 77 > independent of a. 
Combining these estimates, we obtain 



e^^PxdX 



<C \Xry2 ^ReXt^eW/^\.-y\ ^^ 



(157) <Ce-'^'' R-^l\-'^''^''' Rdi^ 



-arg(Ao) 
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Likewise, 

Jr2 Jt2 

(158) Jr2 

<Ce-^"'* [ |/mA|-l/2g-,?|/mA-/mAoIt l^j^^l 

Jr2 

Combining these last two estimates, we have 

(159) /// < C7t-V2e-e°'*/2g-(x-j/)V8et < ^^-i/2g-„ig-(x-y)V8et^ 

for r/ > independent of a. Observing that \x — at\/2t < \x — y\/t < 
2\x — at\/t for any bounded a, for \x — y\/t sufficiently large, we find that 
/// may be absorbed in any summand t~^''^e~^^~y~"'k^> /™*. 

Term IV . Similarly, as in the treatment of the term ///, the principle value 
integral for the "parabolic error term IV may be shifted to r] = R = 9a'^, a 
as above. This yields an estimate 

/ + 00 
-oo 

absorbed in ©(e-^^e-l^-^^l'/^-^*) for all t. 

Case II. \x — y\/t bounded. We now turn to the critical case where 
\x — y\/t < S, for some fixed S. 

Decomposition of the contour: We begin by converting the contour inte- 
gral (I152p into a more convenient form decomposing high, intermediate, and 
low frequency contributions. 

We first observe that L has no spectrum on the portion of 0, lying outside 
the rectangle 

(160) ^e := {A : -r/i < 3f?A < r/, -i? < 9A < R} 

for r] > 0, R > sufficiently large, hence G\ is analytic on this region. 
Since, also, H\ is analytic on the whole complex plane, contours involving 
either one of these contributions may be arbitrarily deformed within 0,\^e 
without affecting the result, by Cauchy's theorem. Likewise, Pa is analytic 
on Qp \ 3f?e, and so contours involving this contribution may be arbitrarily 
deformed within this region. Thus, we obtain 

Observation 3.1 fj MaZ3J ). Assume (H0)-(H3) and (D) hold. The princi- 
ple value integral (|152p may he replaced by 

(161) G(x, t; y) = Ia + h + Ic + Ha + Ih + HI 
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where 






la 


l"rj+ioo 

■.= P.V. e^'Hxix,y)dX 

J Tj—ioo 




h 


■■=P-y{ + je'^'iGx-H^- 


-Px){x,y)dX 


Ic- 


= [ e^*PA(x,2/)dA 




Ila 


, f-rii-ir/2 f-rji+iR , 

:=(/ +/ )e^'Gx{x,y)dX 




lib 


r-rii+iR 

■■=- e^'Hx{x,y)dX 

J —T]i—iR 




III 


:= [ e^'Gx{x,y)dX 





with 
Ti : = [— T/i — ir/2, r] — ir/2] U [f] — ir/2, rj + ir/2] U [77 + ir/2, —rji + ir/2] 

T2 :=dnp \ n, 

for any 7],r > 0, R sufficiently large, and rji sufficiently STnall with respect 
to r. 

Using the above decomposition (J16ip . we shall estimate in turn the high- 
frequency contributions la^h, and Ic, the intermediate-frequency contribu- 
tions Ha and lib, and the low- frequency contributions ///. 

High-frequency contribution. We first carry out the straightforward esti- 
mation of the high-frequency terms la, h, and Ic- The principal term la has 
already been computed in (jl53p to be II{x, t; y). Likewise, calculations sim- 
ilar to those of Term II show that the term lb is time-exponentially small. 
For example, the term e^^X~^B{x,y;X) contributes 

-m{t-i; 1/A4z)dz)^~ f;{^,/A.)iz)dz 






(164) e''' ^y ^M* W'^^g- f;{v*/A,Kz)dz ^ (j^ Jil"A,%) ^-e\x-y\ < (j^-e\x~y\/2^ 

for r/i sufficiently small. This contributes in the term C)(e~''^*^*'^'^^^'') of R. 
Likewise, the contributions of terms e'^*A^^(x— y)C(x, y; A) and e'^*A~^Z)(x, y; A) 



(162) 


X e-'^'^'-Jy ^/^ 


where 




(163) 


P.VA / + / 

^J-00 JR 


and 
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split into the product of a convergent, uniformly bounded integral in ^, a 
bounded factor analogous to (jl64p . and the factor e"''^*, giving the result. 

The term Ic may be estimated exactly as was term III in the large \x—y\/t 
case, to obtain contribution C'(t~^'^e~''^*) absorbable again in the resid- 
ual term C'(e~''*e~l^~^l '*^*) for t > e, any e > 0, and by any summand 

C'(t-l/2(l + t)-l/2e-(x-?/-a+)2/Mt)g-^(^ ^ y) fo^ ^ ^^^11 

Intermediate-frequency contribution. Error term //f, is time-exponentially 
small for rji sufficiently small, by the same calculation as in (J162p - (jl64p . 
hence negligible. Likewise, term Ha by the basic estimate ()64p is seen to be 
time-exponentially small of order 0{—r]it) for any r/i > sufficiently small 
that the associated contour lies in the resolvent set of L. 

Low-frequency contribution. It remains to estimate the low-frequency term 
///, which is of essentially the same form as the low-frequency contribution 
analyzed in |ZH1 lYZJ in the strictly parabolic case, in that the contour is the 
same and the resolvent kernel G\ satisfies same bounds (with no Ex term) 
in this regime. Thus, we may conclude from these previous analyses that 
/// gives contribution as claimed, exactly as in the strictly parabolic case. 
For completeness, we indicate the main features of the argument here. 



Bounded time. For t bounded, we can use the medium-A bounds |Ga|, 
|Ga^|, \Gxy\ < C to obtain | Jp e''^^GxdX\ < C2|ri|. This contribution is 
order Ce~^^ for bounded time, hence can be absorbed. 



Large time. For t large, we must instead estimate Jp e G\dX using the 
small-|A| expansions. First, observe that, all coefficient functions djk{\) are 
uniformly bounded (since |A| is bounded in this case). 

Case II(i). {Q < y < x). By our low- frequency estimates in Proposition 
12.211 we have 

e^'Gx[x,y)dX= f J2e^'d,k(p+{x)i^+{y)dX 
1 -'ri ■ ■ 



(165) 






where each djk is analytic, hence bounded. We estimate separately each of 
the terms 

/ e^'d,,cl)+{x)i^+iy)dX 
on the righthand side of (jl65p . Estimates for terms 

e^'cP+{x)4>t{y)dX 

1 
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go similarly. 



Case Il(ia). First, consider the critical case a^ > 0, at < . For this 



case, 



where 



Set 



4{X) = -X/al + A2/3+/(a+)3 + 0{\^) 
p+(A) = -A/a+ + A2/?+/(a+)3 + 0(A3). 



a = L_ , p.,= _ <o. 

Define F'^^^^ to be the portion contained in VLq of the hyperbola 
(166) 
Re{p+x - v+y) + 0(A3)(|x| + |y|) 

= {l/al)Re[K-alx/a+ + y) + X^xP+at /{a+f - y^t /(<?)] 

= constant 

= (l/a+)[(A™„(-a+x/a+ +y) + Xl,JxP+a+ /{a+f - y/3+/(a+)2)], 

where 

{^ if 1^1 < e 

Denoting by Ai, AJ, the intersections of this hyperbola with dO,g, define 
F'^^ to be the union of AiAq and AqA^, and define F'^^ = F'^^^ UF'^^. Note that 
A = a/p minimizes the left hand side of (I166p for A real. Note also that that 
p is bounded for a sufficiently small, since a < e implies that 

(|a+x/a+| + |y|)/t<2|a+| + 2e 

i.e. (|x| + \y\)/t is controlled by a. 

With these definitions, we readily obtain that 

Re{Xt + pjx - v+y) < -{t/a^){a^ /4p) - r]Im{Xft 

(168) „ „ 

< -a^t/M - r]Im{Xft, 

for A G F'^ (note: here, we have used the crucial fact that a controls 
(|x| + \y\)/t, in bounding the error term 0(A^)(|a;| + \y\)/t arising from 
expansion Likewise, we obtain for any q that 

(169) / \X\'^e^<^'+pt''-''^yUx < Ct-^-§e-"'*/^^ 



for suitably large C, M > (depending on q). Observing that 

a = {a+/a+){x-a+{t-\y/at\))/2t, 
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we find that the contribution of ()169|) can be absorbed in the described 
bounds for t > \y/a^\. At the same time, we find that a > x > for 
t < lu/o-kl^ whence 

a>{x-y- ajt)/Mt + \x\/M, 

for some e > sufficiently small and M > sufficiently large. 
This gives 

g-aVbl < ^-{x-y-a+t)yMt^-r,\x\ 

provided \x\/t > a^, a contribution which can again be absorbed. On the 
other hand, if t < |x/a^|, we can use the dual estimate 

a = (-y-a+(t-|x/a+|))/2t 
>{x-y-a+t)/Mt + \y\/M, 

together with \y\ > \a7t\, to obtain 



a contribution that can likewise be absorbed. 

Case Il(ib). In case a^ < or a^ > 0, terms \(f^\ < Ce^"^'^' and 

l^^l < Ce~''l^l are strictly smaller than those already treated in Case Il(ia), 
so may be absorbed in previous terms. 

Case Il(ii) {0 < x < y). The case < x < y can be treated very similarly 
to the previous one; see |ZHj for details. This completes the proof of Case 
II, and the theorem. 

4. Energy estimates 

4.1. Energy estimate I. We shall require the following energy estimate 
adapted from |MaZ4t IZ2| . Define the nonlinear perturbation variables U = 
(n, v) by 

(171) U{x,t):=U{x,t)-U{x). 

Proposition 4.1. Under the hypotheses of Theorem \1.7\ let Uq G H"^ and 
U = {v,u) be a solution of ([2]) and (I17ip . Suppose that, for < t < T, 
the Wx' norm of the solution U remains hounded by a sufficiently small 
constant C > 0. Then 

(172) \\Uml^<Ce'''\\Uo\\l. + cl\-'^'--^{\\U{T)\\l,+B,ir))dT 
for all < t < T , where the boundary operator Bh is defined in Theorem 
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Proof. Observe that a straightforward calculation shows that \U\h^ ~ | Vl^l^/'- , 

(173) W = W -W :=W{U)-W{U), 

for < r < 4, provided |C/|^y2,oo remains bounded, hence it is sufficient to 
prove a corresponding bound in the special variable W. We first carry out a 
complete proof in the more straightforward case with conditions (A1)-(A3) 
replaced by the following global versions, indicating afterward by a few brief 
remarks the changes needed to carry out the proof in the general case. 

(Al') A{W),A'^,A^^ are symmetric, i° > 6*0 > 0, 

(A2') no eigenvector of A{A^)~^{W) lies in the kernel of B{W), 

(A3') W = (^i) ,^ = (° ~) ,6 > ^ > 0, and G = 0. 

Substituting ()173|) into (JH), we obtain the quasilinear perturbation equa- 
tion 

(174) A'^Wt + AW^ = (BW^)^ + MiW^ + {M2W^)^ + g{W,) - g{W^) 

where A^ := ^"^(T^ + W) is positive definite symmetric, A := A{W + W) is 
symmetric. 

Ml = A{W + W)- A{W) = ( f dA{W + 0W)de]w, 

M. = B(W + W) - BiW) = (» ^^, ^^^,^, » ^^^.^^^^;^ . 

As shown in |MaZ4j . we have bounds 

(175) l^"! < C, |A0| < C\Wt\ < C{\W^\ + \wii\) < CC, 

2 

(176) \d,A'\ + \dlA'\ < CiY, \d'.W\ + \W,\) < C(C + \W,\). 

fc=i 

We have the same bounds for A, B, K, and also due to the form of Mi,M2, 

(177) \Mil\M2\<C{C+\W^\)\W\. 

Note that thanks to Lemma [1.31 we have the bound on the profile: \Wx\ < 
Ce~^l^l, as X ^ -l-cx). 

The following results assert that hyperbolic effects can compensate for de- 
generate viscosity B, as revealed by the existence of a compensating matrix 
K. 

Lemma 4.2 ( |KSh| ). Assuming (Al'), condition (A2') is equivalent to the 
following: 

(Kl) There exists a smooth skew- symmetric m,atrix K{W) such that 

(178) ^e{K{A^)-^A + B){W) > ^2 > 0. 
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Define a by the ODE 

(179) a^ = -sign{A^^)c,\W^\a, a(0) = 1 

where c* > is a large constant to be chosen later. Note that we have 
a^/a = I3x/l3 + 7^/7 and thus 

(180) {a^/a)A^^ < -cJi\W^\ =: -u}{x) 
and 

(181) \a^/a\ < c^\Wa:\ = e^'^uj{x). 

In what follows, we shall use (,) as the Q-weighted L^ inner product 
defined as 

if, 9) = {af,g)L^ 

and ll/lls = Ei=o(^/'^/)^^^ as the norm in weighted H^ space. Note 
that for any symmetric operator S, 

{Sh,f) = -\{{S^ + {a,/a)S)f,f) - ^Sofo-fo- 

Note that in what follows, we shall pay attention to keeping track of c*. 
For constants independent of c*, we simply write them as C. 

4.1.1. Zeroth order "Friedrichs-type" estimate. First employing integration 
by parts yields, and using estimates (|175p . (|176p . and then (jlSOp . we obtain 

-{AW^,W) 

= ^{{A,, + {aJa)A)W, W) + ^AoW{0) ■ W{0) 



1 
1 



< ^{{aJa)A^^w',w') + C7((C + \W^\)\W\ + Loix)\w''\,\W\) + 4° 



<--{oj{x)w\w') + C{C\\w'U + {\W,\w',w')) + C{c,)\\w''\\i + ll 

where 1° denotes the boundary term |AoVF(0)-I^(0), and the term {\Wx\w^ ,w^) 
may be easily absorbed into the first right-hand term, as for c^= sufficiently 
large, 

(182) {\W,\w^,w^) < {c,9i)-\ujix)w^,w^) < -^{u;{x)w^ ,w^). 

Whereas, integration by parts also yields 

{{BWx)x, W) = -{BWx, Wx) - {{ax/a)BWx, W) - B^W^ ' W{0) 

< -GW^i^Wo + C{u;{x)wi^, u) - bowi^{0)w^\0) 

< -eWw'JWl + Cic,)\\w''\\l - bowi'iO)w'\0). 

where we used the fact that BqWx{0)-W{0) = bow!j.^ {0)w^ ^ (0) , and similarly, 
together with the block-diagonal form of B and thus of M2 = B{W + W) — 
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B{W), 

m2w,)x.w) 

= -{M2W,,W,) - {{a,/a)M2W,,W) - M2(0)VF,(0) • VF(0) 

< C(|VF^||VF|, |wf I) + C{lo(x)\W\,w^^) - M2{0)W^{0) • T^(0) 

< ^Ww'JWl + C(^e{u;{x)w',w') + C{c,)\\w''\\l^ - M2{0)W,{0) ■ W{0) 

for any small ^. Note that C is independent of c*. Therefore, for ^ = 9/2 
and c^ sufficiently large, we can compute 

(183) ~{A'W,W) 

= {A''Wt,W) + ^{A',W,W) 

= {-AW^ + {BW,), + M^W, + {M2W,)^,W) + ^{A^tW, W) 

(184) < -^[{u;ix)w',w') + e\\w'J\\l] + CCWWl + C{c,)\\w''\\l + 4° 
where the boundary term 

(185) 4° := ^AoW{0) • W{0) - bow'J {0)w" (0) - M2{0)W^{0) • W{0) 

which, thanks to the negative definite of Aq = A^^{0) and the special form 
of M2 , is estimated as 

(186) 4° < -^-^\w\or + c{\w''{or + \wi'm\nj''m) 

for the outflow case, and similarly 

(187) i^,<c{\wm' + WJm\y^''m) 

for the inflow case. 

Therefore together with this, (jl84p yields 

(188) <A[l_u^x)w',w')+0\\w''\\l] + CX\\w'\\l + C{c,)\\w''\\l + I°. 

4.1.2. First order "Friedrichs-type" estimate. Similarly as above, we shall 
need the following key estimate. We compute by the use of integration by 
parts, (jl82p . plus c^ being sufficiently large, 

-{W^„ AW,,) = ^{W,, {A, + {a,/a)A)W,) + ^AoW,{0) ■ W,{0) 

(189) < -^{oo{x)wi,wi) + CCWwiWl + CcllWJWl 

(190) +^AoW,{0)-W,{0). 
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We shall deal with the boundary term later. Now let us compute 

-1 7 -1 

(191) 2dt^^°^"'^"^ " (T^x,(^°m)x) - {W^,AlWt) + -{A'lW^,W^ 



We shall control each term in turn. First we have 

x\\o 



{A'lW^,W^)<CC\\W^''^ 



and by multiplying {A^) "^ into (I174|) . 

\{w^,AlWt)\ <c{ic + \mw,\,{\w,\ + \wii\ + \W\)) 

mwii\\i + c{{c+\w,\)wi,wi) 

+ C{{( + \W,\)w',w') + C\\w''\\l 

where, once again, the term {\Wx\wl.,w^) may be treated the same way as 
(jl82p . Using ()174p . we write the first right-hand term in ()19ip as 

{Wx, {A'>Wt)x) ={W„ [-AW, + {BW,)x + M^W, + (M^W^U^) 
= - {Wx, AWxx) + {Wx, -A^Wx + {M^Wx\) 

- (t^xx + {axloL)Wx, \{BW^\ + (M2l^,.)x]) 

- WM\iBWx\ + (M2H^JJ(0) 
1 



< 

- 4L 

+ C7 






^ll 



where l\ denotes the boundary terms 



(192) 



1 



-A^WM ■ Wx(0) - WM ■ \{BWx)x + (M2Pr.).](0), 



and we have used estimate (|190p . (|182p for w^jwj^, and applied Young's in- 
equality to obtain: 



{Wx, -A^Wx + (MiW^).,) < cue + \W,\)\Wx\,\W,\ + \W\). 

-0\Ui\\l 

+ C{\wii\+u{x)\wi'\,{C + \Wx\)\wi'\) 



-{Wxx + {ax/a)Wx, {BW^).,) < -9\\w'JX 



< 



-{Wxx + {ax/a)Wx,{M2Wx 

C{\w'J,\ + ^(x)l^f I, (C + \Wx\){\Wx\ + \W\)). 
Putting these estimates together into (|19ip . we have obtained 



(193) 



< C 



II\\2 



CW\\i + {\W,\w\w') + C{c.)W\\i 



+Il 
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Let us now treat the boundary term. First observe that using the paraboUc 
equations with noting that A^ is the diagonal-block matrix diag(A]') ^2)) '^^ 
can write 

u;f (0) • [A^wl' + A2iwi + A22wi^ - MiI^J(O) 

< e\wi'm' + c{\wm' + \w'jm' + i^f'wn. 

Whereas, for the first term in I^, we shall consider two cases separately. 
When A^^ < -Oi < 0, we get 



AoW^iO) ■ WM < -y l^^(0)P + C\wi'm'- 



Therefore 



(194) II < --^\wim' + C{\W{Or + kf (0)P + \wi'm')- 

Meanwhile, for the case A^^ > ^1 > 0, we have 

\AoW^{0)-W^{0)\ <C\W^{0)\^. 

Now the invertibility of A^^ allows us to use the hyperbolic equation to 
derive 

\wim<ci\wim + \wi'm)- 

Therefore in the case of A^^ > Oi > 0, we get 

(195) /I < c{\wm' + mm' + i^f (o)p). 

Now applying the standard Sobolev inequality (applies for a— weighted 
norms as long as |ax/a| is uniformly bounded): 

(196) \w{0)\' < C\\w\\L2i\\w,\\L2 + 1^11^2) 
to control the term |tt;^^(0)p in I^ in both cases. We get 

(197) WJ{0)f <e'\\wii\\l + C\\wi'\\l 

Using this with e' = 6/8, (fT92l) . and (fTM]) . the estimate (fT93]l reads 

(198) <c(^cW\\l + {\W,\w',w')+C{c,)\\w''\\j)+ll 
where the boundary term I^ is estimated as 

(199) II < -|k^(0)|2 + C{\W{Or + \wi''-'^'' 
for the outflow case, and similarly 

(200) ii < ci\wm' + mm') 
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for the inflow case. 



4.1.3. Higher order "Friedrichs-type" estimate. Similarly as above, we shall 
now derive an estimate for {A'^d^W, d!^W), k = 2,3, 4. We need the following 
key estimate. Integration by parts and (I180p give 

-{d',W,Ad^,~^'W) =^(5>, {A, + {a,/a)A)d^,W) + ^Aod'^WiO) ■ d'^WiO) 

We compute 

^±{A'd',W,d',W) =^{A',d',W,d',W) + {A''d',W,d',Wt) 

=^{A',d',W,d',W) + {A''d',W,d'JiA''r' 

(201) i-AW, + (BW,),) + MiW^ + {M2W,),]). 

We shall estimate each term in turn. First, \{A'^dl^W,d^W)\ < CC\\d^W\\l, 
and 

{A^d^W,d^[-{A^)~^AW^]) 

k 

= {A^d',W,^di,[-{A^)-'A]dt'^'W) 

k 

= -{d^w,Ad^+^w) + ^{A°d^w,di[-{A^)-^A]dt'^^w) 

i=\ 

where we have 

(202) |5;[-(A0)-1A]|<C7 Y^ n l^-^'^l- 

Using the hypothesis on the boundedness of solutions in W'^'°^ and weak 
Moser inequality [Z4i Lemma 1.5], we get 

\{A'^dlW,cf,\-{A'^)-^A\dl-^^^W)\ < 

k 

c[\\w''\\l + C\\w'\\l + Y{\w,\dW,dW))- 

i=l 
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This, dSnU), similar treatment ([182]) for {\Wa:\dl;w^ ,dl;w^) with c* being 
sufficiently large give 

{A''d',W,d'A-{A')-'AW,]) < -i(a;aV,5.'^^) + ^^o^'W^W • d',W{0) 

fc-i 

(203) + c[\\w''\\i + cww'wi + J2{mdw, dW] 

i=l 

Next, similarly, we obtain 



Finally, we compute and 
{A''d',W,d',[{A'^)-\BW, + M2W,U) 

k 

= ^{A^d^W,di[{A°)-^]dt'^\BW^ + M2W^)) 

i=0 

={d^W,d^+\BW, + M2W,)) 

k 

+ Y,{A''dlW,di[{A'')-^]dl-'+\BW^ + M2W^)) 

i=l 

< - {d^+^W + ia,/a)d^W, d^{BW, + M2W,)) 

k 

+ ^\\d'.^'u,''\\l + C[cl\\w''\\l + CWw'Wl + Y,{\W.\dW,diw')) 

i=l 
k 

+ C{cl\\w^'\\l + CWw^Wl + ^{mdlw^,dlw^) 

where in the last inequality we used the special form of B and M2 to get 
{d^+^W + {a^/a)d^W,d^{BW^, + M2W,)) 

k 

< -9\\d',+W\\l + C[c{c.)\\w''\\l + Clk^ll^ + Y^imoW, dW) 

i=l 

Note that in the last inequality, there is no term of {uj{x)dl.w^ , d^w^) because 
of the presence of \Wx\ in term of 112 M2. 



STABILITY OF BOUNDARY LAYERS 43 

Put all these estimates into (120 ip together, we have obtained 
l±{A^dlW,dtW) + \e\\dl+W\\l + \{u:{x)dW,dW) 

k~l 

(204) <c(c{c.)W\\l+C\W\\l + Y.{\W,\dW:dW))+h 

i=l 

where the boundary term 

(205) h := -AodlW{Q) ■ a^W(O) - d^^[hd.^w^^ + M|2vK,](o)5^u;^^(0). 

For this boundary term, we shall treat the same as we did before. First 
using the parabolic equations with noting that AP is the diagonal-block 
matrix diag(j4']*,^2)5 ^^ ^^"^ write 

(206) =d^,-^\Al{Q)wi\Q,t)+A2iw'M+A22w'J{Qi)-Il2Mi{Qi)WM]. 
Therefore we get 

k 

< c\dW'm [\d^.-'w{'m + ^(ia>^^(o)i + \diw\m 

1=0 

k k 

2 



(207) <eY,\dWm'+cY,\dWm 

i=0 i=l 

(208) +C\d',w''m\dtWm 



for any e small. To deal with the term of wj.^ , for simplicity, assume k = 3. 
By solving the parabolic-part equations and using the invertibility of b, we 
obtain 

(209) ., 

\W^t\ < C{\W\ + \w,\ + \w,,\ + \wiU). 

Since for case A: = 3 we have a good term ||9^i(;-'^-'^||o (see ()204p ). the 
term |ty^;^a;(0)| resulting from the boundary treatment is easily treated via 
Sobolev embedding inequality. Hence all terms in a form dJ^w (0) are easily 
estimated. Meanwhile, using the hyperbolic-part equations, we have 

(210) \w(\<C{\W\ + \W^\). 

Employing Young's inequality to the last term in (j207p . we obtain 

(211) < e ^ \dW{o)f + ciY^ \diw''m' + \wi'm' + K'mi') 

i=0 i=0 
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To deal with the term of w^ , we need to consider two cases separately. When 
^" < -9i < 0, we get 



Therefore 



Aod'^WiO) ■ d'^WiO) < -^|9V(0)|' + C|aV 



k—l 

It < - y |5V(0)|' + C(^ \diw\0)\' 

k 

(212) + ^ ia>^^(o)p + \wi'm'' + \wiim')- 

1=0 

Meanwhile, for the case A^^ > ^i > 0, we have 

I^O^XO) • 5X0)| < C(|5V(0)P + |5V(0)P). 
The invertibility of ^^^ allows us to use the hyperbolic equation to derive 

i=o 
Therefore in the case of A^^ > 9i > 0, we get 

k 

(213) It < c{j2{\dWm' + \diw'm') + \wi'm' + \wiim 



Employing the boundary estimates into (j204p . we have obtained 



k—l 

(214) <c[c\\w'\\l + cl\\w''\\l + J2{\WMy^':diw')) +It 

j=0 

where, after absorbing the terms of \dj^w^^{0)\ via Sobolev embedding 
boundary term I^ satisfies 

k—l 

(215) It < -^-^idtw'm' + c(^ \dWm\' + kr(o)p + \w{!m' 

1=0 

for outflow case, and 

k 

(216) It < C(^ \diw\0)\' + \wi\0)\' + \wilm') 

for the inflow case. 

\k 



We shall establish an Kawashima-type estimate to bound the term ||u;^"^ 



appearing on the left hand side of the above. 
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4.1.4. "Kawashima-type" estimate. Let K be the skew-symmetry in (I178p . 
Integration by parts and skew-symmetry property of K yield 

{KW^t, W) = -{KWt, W^) - {{K,, + {a^/a)K)Wt, W) - KqWo • {Wo)t 
= {KW,., Wt) + {{K, + {aJa)K)W, Wt) - KqWo ■ iWo)t. 

Using this, we compute 

{KtW^ + KW^t, W) + {KW^, Wt) 
= {KtW,, W) + {2KW, + {K, + {a,/a)K)W, Wt) 

- KoWo ■ {Wo)t 
= {KtW,, W) + {2KW^ + {K, + {a^/a)K)W, -{A^)-^AW^) 

+ {2KW^, + {K, + {a,/a)K)W, {A'')-\BW^), 

+ M,W^ + {M2W^)^) - KoWo ■ {Wo)t 
<-2{K{A'^)'^AW^,W,)+C\\wi\\l - KoWo ■ {Wo)t 

+ C(C{C,)\W'\\1 + Q\\w'\\l + {u:{x)w' y) + {u{x)wlwi)). 

Using ()178l) . we get 

{K{A'')-^AW,,W,) > 92\\wi\\l - Cico)\\wi'\\l 
and thus obtain from the above estimate with ^ = 62/2 

j^{KW,.W)<-^\\wi\\l + c(c{c.)\W'\\l + C\W\\l 
(217) + {oj{x)w',w') + {uj{x)wlwi)) - KoWo ■ {Wo)t. 



4.1.5. Higher order "Kawashima-type" estimate. With similar calculations, 



we shah obtain an estimate for 4f{Kd^W.,d^ ^W).,k > 1. We compute 



{Kd^Wt,d^-^W) = {Kd^W,d^-^Wt) 

+ {{K, + {a,/a)K)dt'W, d^'Wt) - Kd^'Wt ■ d^'WiO). 
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and hence 

j{KdlW,dl-^W) = {Ktd^^W,dl~^W) + {2Kd';W,d';-^Wt) 

+ ((i^, + {a^/a)K)dl-^W,d^-^Wt) - Kd^'Wt ■ dt'W{0) 

+ {{K^ + {a,/a)K)dt^W, 

- Kd^-^Wt ■ d^-^WiO) 
< -2{K{A''r'Ad',W, d'^W) + e\\w'\\l + Ccl\\w''\\l_,, 

k 

+ CCWw'Wl + C^{uix)diw\ dW) - Kdt'Wt ■ dt'W{0) 
1=1 

for e small. 

Using (|178l) . we obtain from the above 

(218) §{Kd'.W, dt'W) < - |||9Vllo + Ccl\\w''\\l^, + 6||«;^||,_i 



+ CC\\w'\\l + CY,Hx)dW,dW) 
1=1 

(219) -Kd^-^Wfd^-^W{0). 



4.1.6. Final estimates. We are ready to conclude our result. First combining 
the estimate (jl98p with (|188p . we easily obtain 

i| ((A°W,, VT,) + M(A0VF, t^)) 

<-{^\\wii\\l + ^{^{x)wi,wi)^ 

+ C[c\\w'\\l + {mw',w') + C{c,)\\w''\\l) + II 

-^(^{u;{x)w',w')+9\\wi'\\l)+CMC\\w'\\l 

+ MC{c,)\\w^^\\l + MI^ 

By choosing M sufficiently large such that M9 ^ C{c^), and noting that 
c*0i|Wa;| < ^{x)i we get 

^^^(^{A'^W.^W,) + M{A''W,W)) 

(220) < - [e\\w^^\\l + {u{x)w^, w^) + {uj{x)wi,wi)'^ 

+ c[c\W\\l + C{c.)\\w''\\l) + 4^ + Mil 
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We shall treat the boundary terms later. Now we employ the estimate (I217p 
to absorb the term ||w ||i into the left hand side. Indeed, fixing c^, large as 
above, adding (|220p with (j217p times e, and choosing e, C, sufficiently small 
such that eC{c^,) ^ 0, e ^ 1 and C, <^ €02, we obtain 

i^ ((^°H^,., W^) + M{A^W, W) + e{KW^„ W) 



< - (eww^^wl + {uj{x)w^,w^) + {uj{x)wi,wi)) 



/||2 I r^f„ \ii„.JIii2\ ^se,, , r||2 



'"l + W'\\l)+h 



+ Ce(c(c,)||t/;^^||i + CWw'Wl + (a;(x)t/;^^') + {u;{x)wi,wi) 

+ /I + M/O - eKoWo ■ {Wo)t 

<-l{0Wf2 + e2e\\wi\\l)+C{c,)[C\\w 

where h := 4' + M/O - eKoWo ■ iWo)t. 

By a view of boundary terms 1)1,11, ^^ treat the term /^ in each in- 
flow/outflow case separately. Recalling the inequality (jl97p . |i(;^^(0)| < 
C||w II2. Thus, using this, for the inflow case we have 

h < M\wm^ + cim(o)|2 + M\wi'm\w''m 

< ^\\w''\\l + M^\W{0)\^ + C\Wt{0)\^. 

Meanwhile, for the outflow case, with M61 » 1 and KQWQ-{WQ)t ~ Wq^Wq^A- 
WqWqI, we have /& is bounded by 

+ C{\wi\G)\'' + k^^(O)p) + 6(|u;f (0)p + Wtm-") 
which, together with e being sufficiently small and the facts that 

\wim<c{\wim + WJm + \wm) 

obtained from solving the hyperbolic equation and the embedding inequality 

\wi\Q)\<C\Wh, 
yields 

h < -|(ki(o)p + k^(o)p) + ^-W^wl + c{\w^^{or + \w{\0)\') 

for the outflow case. Now by Cauchy-Schwarz's inequality and by positivity 
definite of A^, it is easy to see that 

(221) £ := {A°W^,Wa:) + M{A^W,W) +e{KW^,W) r. WWfrr^ -- \\Wfrrr. 

The last equivalence is due to the fact that j3 is bounded above and below 
away from zero. Thus the above derives 
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j^£{W){t)<-9s£iW){t) + Cic.){\\wml, 

for some positive constant ^3, which by the Gronwall inequaUty yields 
(222) 



WWml, < Ce-''\\Wo\\j,.+Cic.)|\''('-^\\\W{T)\\l,+B^{T))dT, 



where W{x,0) = Wo{x) and 

(223) B^{T) := O{\W{0,t)\^ + |i^t(0,r)p) = 0{\{h,,h2)\^ + |(/ii, /i2)tP) 
for the inflow case, and 

(224) ^i(r) := 0{\w'' {0,t)\^ + \wi'{0,T)\') = 0{\h\' + \ht\') 

for the outflow case. 

Similarly, by induction, we shall derive the same estimates for W in H^. 
To do that, let us define 

£i{W) := {A^W^, W^) + M{A°W, W) + e{KW^, W) 
£k{W) := {A°d^W,d^W) + M£k-i{W) + e{Kdl;W,d^~^W). 

Then by Cauchy-Schwarz inequality, it is easy to see that £kiW) ~ 
llVFll^fc, and by induction, we obtain 

j/kiwm < -e3£k{w)it) + c{c,)i\\wit)\\l, + Bkit)), 

for some positive constant ^3, which by the Gronwall inequality yields 

(225) 

Jo 
where W{x, 0) = Wo{x), k = 2,3, 4, and Bk are defined as in (fTOll and (1201) . 

4.1.7. The general case. Following |MaZ4j . the general case that hypothe- 
ses (A1)-(A3) hold can easily be covered via following simple observations. 
First, we may express matrix A in (|174p as 

(226) AiW + W) = A + {C + m\)(^Jl^^ ^[J] 

where ^ is a symmetric matrix obeying the same derivative bounds as de- 
scribed for A, identical to A in the 11 block and obtained in other blocks jk 
by 

A^\W + W)= A^''{W) + A^'^iW + W)- A^\W) 

(227) = A^'^iW^) + 0{\W^\ + \W^\) = A^^{W+) + 0{C + \W,\). 
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Replacing ^ by A in the /c*^ order Priedrichs-type bounds above, we find 
that the resulting error terms may be expressed as 

plus lower order terms, easily absorbed using Young's inequality, and bound- 
ary terms 

k 

o(£\d>"m\dWm) 

resulting from the use of integration by parts as we deal with the 12— block. 
However these boundary terms were already treated somewhere as before 
(see (I207p ). Hence we can recover the same Friedrichs-type estimates ob- 
tained above. Thus we may relax (^1') to (^41). 

The second observation is that, because of the favorable terms 

occurring in the lefthand sides of the Friedrichs-type estimates ()214p . we 
need the Kawashima-type bound only to control the contribution to \d^w {"^ 
coming from x near -|-oo; more precisely, we require from this estimate only 
a favorable term 

-02((i-o(c + |i^.l))aV,9V) 

rather than 02||c^^w^|lo ^s in (j218p . But, this may easily be obtained by 
substituting for K a skew-symmetric matrix- valued function K := K{W+), 
and using the fact that 

^e{K{A^)-^A + B){W+) > ^2 > 0, 

and same as (|227ll . K = K + 0{C + \W^\), we have 

^eiKiA'^y^A + B){W) > ^2(1 - 0(C + I^^^D) > 0. 

Thus we may relax {A2') to {A2). 

Finally, notice that the term g{Wx) — g{Wx) in the perturbation equation 
may be Taylor expanded as 

\ f ^ 

gi{Wx,Wx)+gi{Wx,Wx)) + [oi\Wx\' 



The first, since linear term on the righthand side decays at plus spatial 
infinity and vanishes in the 1-1 block, it may be treated as follows by Young's 
inequality 

\ ^''}^ <c({{C + \Wx\)wi,wi) + \\wi'\\l 



giiWx,Wx)+giiWx,Wx)J \w^ 

which can be treated in the Friedrichs-type estimates. The (0, 0(|Wa;p) 
nonlinear term may be treated as other source terms in the energy estimates. 
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Specifically, the worst-case term 
may be bounded by 

\\dl-^W'\\LA\w\\w^M\w\\H^-dW\^)dl-^o{\w,\^m. 

The boundary term will be contributed into the form (j205p of /fe, and hence 
using the parabolic equations to get rid of this term as treating in (|206p . 
Thus, we may relax (^43') to (^43), completing the proof of the general case 
[Al) — (j43) and the proposition. D 

4.2. Energy estimate II. We require also the following estimate: 

Lemma 4.3 ( [HR] ) . Under the hypotheses of Theorem \l. 71 let Eq := ||(1 + 
|xp)^'^C/o||j:/4, and suppose that, for < t < T, the W'^'°° norm of the 
solution U of (I23ip remains bounded by some constant C > 0. Then, for 
allO<t<T, 

(228) 11(1 + \x\y/^Uix,t)\\j,, < MSoe^*. 

Proof. This follows by standard Priedrichs symmetrizer estimates carried 
out in the weighted H'^ norm. D 

Remark 4.4. An immediate consequence of Lemma \4-3[ by Sobolev em- 
bedding: W^'°° C -ff^ and equation (1231 p . is that, if Eq and \\U\\}{4 are 
uniformly bounded on [0,T], then 

(229) (1 + |x|)3/2 ^\U\ + \Ut\ + \U,\ + \U,t\\ (x,t) 
is uniformly bounded on [0, T] as well. 

5. Stability analysis 

In this section, we shall prove Theorems 11.61 and 11.71 Following [HZt 
IMaZ3] ■ define the nonlinear perturbation U = {u,v) by 

(230) U{x, t) := U{x, t) - U{x), 
we obtain 

(231) Ut-LU = Q{U,U,),, 
where linearized operator 

(232) LU := -{AU), + (BU,)^ 

where 

AU := dF{U)U - {dB{U)U)U^, B = B{U) 
and the second-order Taylor remainder: 

Q{U, U^) = F{U + U)- F{U) + A{U)U + {B{U + U)- B{U))U^ 
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satisfying 



\Q{U,U^)\<C{\U\\U^\ + \U\'') 

(233) \^^QiU,U^U <Cm\U.\ + \U?) 

\Q{U.U^U<C{\U\\U^^\ + \U,\'' + \U\\U^\) 

so long as |C/| remains bounded. 

For boundary conditions written in f7— coordinates, (B) gives 

h{t) = Ht) -h = iW{U + U)- WiU)){0, t) 

= {dW/dU)iUo)U{0,t) + O{\U{0,t)\'). 

in inflow case and 

h{t) = h{t) -h= {w^\U + U)- w^^{U)){0, t) 

= {dw''/dU){Uo)U{0,t)+0{\U{0,t)\^) 

= m{bi 62)(C7o)C/(0,t) + 0(|C/(0,t)|2) 

= mB{Uo)UiO,t) + Oi\UiO,t)\^). 

5.1. Integral formulation. We obtain the following: 
Lemma 5.1 (Integral formulation). We have 

POD 

U{x,t)= G{x,t;y)Uoiy)dy 
Jo 

+ / (Gy{x, t - s; 0)BU{0, s) + G{x, t - s; 0)AU{0, s)) ds 
(236) 

H{x, t - s] y)UiQ{U, Uy)y{y, s) dy ds 



(235) ,, r ...-....,. .. .....,„ ..,2. 



Jo 

rt /»oo 

- Gy{x,t- s;y)U2QiU,Uy){y,s)dyds 

Jo Jo 

where U{y,0) = Uo{y). 

Proof. From the duality (see |ZHt Lemma 4.3]), we find that G{x,t — s;y) 
considered as a function of y, s satisfies the adjoint equation 

(237) {ds-Ly)*G*{x,t-s;y) = 0, 

or 

(238) -Gs- iGA)y + GAy = {GyB)y. 

in the distributional sense, for all x,y,t > s > 0, where the adjoint operator 
of Ly is defined by 

(239) L;v:=v;A + iv;B)y, 

with V* = V^\ 
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Likewise, for boundary conditions, we have, by duality 

(iii') for all x,t > 0, G{x,t;0) = in the outflow case A^< < 0; and 
G{x, t; 0)B = in the inflow case ^* > 0, noting that no boundary con- 
dition need be applied on the hyperbolic part for the adjoint equations in 
the inflow case. 

Thus, integrating G against (j23ip . we obtain for any classical solution 
that 

ft />oo 

G{x, t - s; y)Q{U, Uy)y{y, s) dy ds = 

'0 ^0 

(240) /■* /"°° 

/ / G{x,t- s;y){ds- Ly)U{y,s)dyds 
Jo Jo 

= :h+h. 

Integrating by parts and using the boundary conditions (iii') on the 
boundary y = 0, we get 

pt /'OO 

h= I / G{x,t- s;y)dsU{y,s)dyds 
Jo Jo 

i-t />oo 

dsG{x, t- s; y)U{y, s) dy ds 

10 JO 

+ G{x,0;y)U{y,t)dy- Gix,t-y)U{y,0)dy 

Jo Jo 

where note that 

Uix,t)= / Gix,0;y)Uiy,t)dy 
Jo 
and also 

pt poo 

h= I / G{x,t- s;y){-Ly)U{y,s)dyds 
Jo Jo 

ft fOO 

G{x, t - s; y)iiAU)y - {BUy)y){y, s) dy ds 
'0 Jo 

ft roo 

-GyA - {GyB)y)U{y, s) dy ds 
/o ^0 

- / Gy{x,t-s;0)BU{0,s)ds- / G{x,t - s;0)AU{0,s)ds 
Jo Jo 

Combining these estimates, and noting that GyB = GyB since HB = 0, 
we obtain ()236p by rearranging and integrating by parts the last term of 

ft roo 

G(x, t - s; y)Q{U, Uy)y{y, s) dy ds 

t rOO 

/ (H + G){x, t - s; y)Q{U, Uy)y{y, s) dy ds 
Jo 

D 



(241) -^o -^o 
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As an expression for Ux, we obtain the following. 

Lemma 5.2 (Integral formulation for Ux)- We have 
(242) 

Uxix,t)= Gxix,t;y)Uo{y)dy- H{x,t - s;0)UiQ{U,Uy)yiO,s) ds 
Jo Jo 

+ / \Gxy{x, t - s; 0)BU{0, s) + Gx{x, t - s; 0)AU{0, s)] ds 
Jo L J 

+ / / {Hx-Hy){x,t-s;y)UiQ{U,Uy)y{y,s)dyds 
Jo Jo 

ft f'OO 

- H{x,t- s;y)UiQ{U,Uy)yy{y,s)dyds 

Jo Jo 

t — 1 POO 

/ Gxyix,t- s;y)U2Q{U,Uy){y,s)dyds 
Jo 

t />oo 

/ Gx{x,t- s;y)Il2Q{U,Uy)y{y,s)dyds 
it-i Jo 

where U{y,0) = Uo{y). 

Proof. Differentiating the formulation (j236p for U{x, t) with respect to x 
and noting that 

i-t roo rt roo 



rt roo rt roo 

/ / Hxc^dyds= / / {Hx-Hy)c^dyds 
Jo Jo Jo Jo 



t roo rt 

j H{x,t — s;y)(f>y{y,s) dy ds — / H{x,t — s;0)(l){0,s)ds 
Jo Jo 



and 



t roo f'^^^ f'°° 

/ Gxyi^ dyds = / Gxyi^ dy ds 

Jo Jo Jo 

t roo rt 

/ Gxipydy ds — / Gx{x,t — s;0)ip{0,s)ds 
t-i ^0 Jt-i 

are valid for any smooth functions (j), tp, we obtain the lemma. D 

5.2. Convolution estimates. To establish stability, we use the following 
lemmas proved in [HZ. HR , RZ] . 

Lemma 5.3 (Linear estimates I). Under the assumptions of Theorem \1.7\ 

\G{x,t;y)\{l + \y\y'/Uy<C{e + i^i+tP2){x,t), 



(243) '' 

/ \Gx{x,t;y)\{l + \y\)-^/^dy<C{e + 2Pi+^P2){x,t), 
Jo 

and so the latter is dominated by ipi + 1^2, for < t < +00, some C > 0. 
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Lemma 5.4 (Linear estimates II). Under the assumptions of Theore7n \1.7\ 
if \Uo{x)\ + \d^Uo{x)\ < Eo{l + \x\)-^/'^,Eo > 0, then, for some 6 > 0, 

r+co 

/ H{x,t;y)Uo{y)dy < CEoe-'\l + \x\)-^/\ 

(244) -^o 

/ H,{x,t;y)Uoiy) dy < CEoe-'\l + \x\)-''/\ 
Jo 

and so both are dominated by C£'o(V'i + 1P2), for < t < +00, some C > 0. 
Lemma 5.5 (Nonlinear estimates I). Under the assumptions of Theorem 

/•t f + CO 

/ / \Gy{x,t-s]y)\^{y,s)dyds<C{e + 7pi+7p2){x,t), 

(245) /° \ 

^ ' r r+00 

/ / \G^y{x,t-s]y)\^{y,s)dyds<C{0 + 7pi+Tp2){x,t), 
Jo Jo 

for < t < +00, some C > 0, where 

(246) *(y,s):=(0 + V'i+V'2)'(y,s). 

Lemma 5.6 (Nonlinear estimates II). Under the assumptions of Theorem 

/ / H{x,t-s;y)T{y,s)dyds<Cii;i+i^2)ix,t) 
Jo Jo 

rt r+OD 

(247) / / iH,-Hy)ix,t-s;y)Tiy,s)dyds<Cii;i+^lj2)ix,t) 
Jo Jo 

30 

\Gx{x,t - s;y)\T{y, s) dyds < C{ipi +ip2)ix,t) 



it~i Jo 
for all < t < +00, some C > 0, where 

(248) T{y, s) := s-'/^O + V'l + V'2)(?/, s) 

We require also the following estimate accounting boundary effects. 
Lemma 5.7 (Boundary estimates I). Under the assumptions of Theorem 

ir2\if\hit)\ + \h'{t)\<Eo{i + t)-\ 

/ H{x, t - s] 0)h{s) ds < CEoiipi + Tp2){x, t) 

(249) ^° 

/ H^{x, t - s; 0)h{s) ds < C7So(V'i + V'2)(x, t), 
Jo 

for < t < +CXD, some C > 0. 
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Proof. Note that H{x, t;0) = for the outflow case A^< < 0. Consider the 
inflow case A^ > (and thus a* > 0). We have 



/ H{x,t-s;0)h{s)ds 
Jo 



1 



a* 



< e-^'ol^lCl + \x- a,t\r^ < CEq{iPi + V2)(x, t), 

/ Hx{x,t — s;0)h{s)ds 
Jo 

<e-''o"/"*f|/i| + |/i'|l(- — (x-a,t))| 
L J a* 

< e-^ol^l(l + |x - aj\)-^ < CEoiiJi + V2)(x,f), 
which completes the proof of the lemma. D 
Lemma 5.8 (Boundary estimates II). Under the assumptions of Theorem 

\r^ if\h{t)\ + \ht{t)\<EQ{i + t)-\ 

'dyix, t - s; 0)Bh{s)+G{x, t - s; 0)Ah{s)) ds 



(250) 



<CEo{e + ^Pl + ^lJ2)ix,t) 



G^yix, t- s; 0)Bh{s)+Gx{x, t - s; 0)Ah{s) ) ds 



<CEoi9 + i^i + i^2){x,t) 



for < t < +CXD, som,e C > 0. 



Proof. The estimate on J^ , where Gy{x,t — s;0),Gxy{x,t — s;0) is non- 
singular, follows readily by estimates similar to but somewhat simpler than 
those of Lemma (j5.5p . which we therefore omit. 

To bound the singular part J^_i, we integrate (I238P in y from to +oo 
to obtain 



(251) GyB + GA 



hoo r+oo 

G{x,t- s;y)Aydy+ Gs{x,t - s;y) dy. 

/o Jo 

Substituting in the lefthand side of (j250p . and integrating by parts in s, we 

obtain 

(252) 



{GyB + GA)h{s) ds 



t-i 



1 , r+oo 

/ Ay{y)G{x,T;y) dy)h{t - t) dT 
^Jo 

1 , /•+00 



+ 



G{x, t; y) dy \h'{t — r) dr 
n'^ G{x,l;y)dy)h{t-l), 



^ JO 

OD 
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which by / \G\dy < C is bounded by maxo<r<i(|/i| + |^'|)(i — t)- 

Combining this with the following more straightforward estimate (for 
large x, |x| > a^t) 
(253) 

/ Gy{x,t-s;Q)Bh{s)ds < I \Gy{x,T;0)\Bh{t - t) dr 
Jt-i Jo 

<C max \h{t-T)\ / T-^e~\^\''/^^ dr 
0<^<l Jo 

= C\x\~'^ max \h{t-T)\ 

0<T<1 

rl 

X / (|x|VT)e-l^l'/^^dr 
< C max |/^(t-r)||x^^ 

0<r<l 

(254) ^ 

/ G{x,t-s]{))Ah{s)ds < / \G{x,T;Q)\Ah{t-T)dT 
Jt-i Jo 

<C max \h{t-T)\ [ r-^/^g-kp/Cr ^^ 

0<r<l' "Jo 

<C max \h{t-T)\ / r-^g-l^l'/^^ dr 

0<T<1 Jq 

< C max |/^(t-T)||x^^ 

0<T<1 

and the estimate (|249p for H term (thus together with (|254p ioi G = G+H), 
we find that the contribution from j _. has norm bounded by 

max {\h\ + \h'\)(t - r)(l + \x\)~^ < CEoii^i + V2)(x,t). 

0<T<1 

Combining this estimate with the one for L , we obtain the first inequality 
in (j250p . For second inequality, we first differentiate ()252p with respect to x 
to get 
(255) 

/ {G^yB + G^A)h{s)ds= U Ay{y)G,,{x,T-y)dy^h{t-T)dT 

Gx{x, t; y) dyj h'{t - r) dr 
(/' G.,{x,l-y)dy)h{t-l), 



'0 ^ JO 

f + OO 

/O 



which, by /q / IG^^Idyfir < G Jq t ^'"^dr < C, is bounded by maxo<r<i(|/i| + 
\h'\){t — t), similarly as above. 
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For the large x, clearly we still have similar estimates as ()253p and (I254p 
for Gxy and Gx- These, estimate ()249p for H^, and ()255p yield the contribu- 
tion from J^_^ as above, which together with the estimate for J„ completes 
theproof of (I250D. D 

5.3. Linearized stability. In this subsection, we shall give the proof of 
Theorem 11.61 We first need the following estimates: 

Lemma 5.9 ( |MaZ4j ). Under the assumptions of Theorem \1.(A 

G{;t-y)fiy)dy < C(l + t)-5(i-V'-)|/|^,, 

/ Hi;t;y)f{y)dy <Ce-^V\Lv, 
Jo LP 

for all t > 0, some G, r] > 0, for any 1 < q < p and f £ L'^ D L^, where 
l/r + l/q = l + l/p. 

Lemma 5.10. Under the assumptions of Theorem \1.6l if \h{t)\ < Eq{1 + 



(256) 







(257) 



/ (Cyix, t - s; 0)Bh{s)+G{x, t - s; 0)Ah{s)) ds 

<CSo(l + t)"5(i-i/p) 
for < i < +00, some C > 0. 



LP 



Proof. This follows at once by the boundary estimate (j250p and the fact 
that 1(0 + ^1 +^2)(-,t)|LP <C(l + t)-5(i-i/p). □ 

Proof of Theorem \1.(A Sufficiency of (D) for linearized stability (the main 
point here) follows easily by applying the above lemmas to the following 
representation for solution f/(x, t) of the linearized equations (jlOp 

/>oo 

U{x,t)= Gix,t;y)Uoiy)dy 
Jo 

+ / (Gy{x, t - s; 0)BU{0, s) + G{x, t - s; 0)^^7(0, s)) ds 

where U{y,0) = Uo{y) and \U{0,s)\ < G\h{s)\ < C(l + s)-^ by (HI]) in the 
inflow case, and \BUiO,s)\ < C|/i(s)| < C(l + s)"^ by ([l2]) in the outflow 
case, noting that G{x, i; 0) = in this case. Necessity follows by a much 
simpler argument, restricting x, y to a bounded set and letting t ^^ oo, 
noting that G is given by the ODE evolution of the spectral projection 
onto the finite set of zeros of D in 3f?A > 0, necessarily nondecaying, plus 
an 0(e~''*) error, rj > 0, from which we find that asymptotic decay implies 
nonexistence of any such zeros; see Proposition 7.7 and Corollary 7.8, |MaZ3j 
for details. D 
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5.4. Nonlinear argument. In this subsection, we shall give the proof of 
Theorem 11.71 In fact, with the above preparations, the proof of nonlinear 
stability is also straightforward. 

Lemma 5.11 {H^ local theory). Under the hypotheses of Theorem \1.7l 
then, forT sufficiently small depending on the H'^—norm oJUq, there exists 
a unique solution U{x,t) G L^{0,T;H^{x)) of (I23ip satisfying 

(258) \U{t)\m < C\Uo\m 

for allO<t<T. 



Proof. Short-time existence, uniqueness, and stability are described in [Z2l 
IZ4] . using a standard (bounded high norm, contractive low norm) contrac- 
tion mapping argument. We omit the details. D 

Lemma 5.12. Under the hypotheses of Theorem\r^ let U G L°°(0, T; H'^ix)) 
satisfy (123ip on [0,T], and define 

(259) C{t):= sup \{\U\ + \U,\){9 + ^Pi + %b2)-Hx,t) 

x,0<s<t L 

If C,{T) and \UQ\m are hounded by (^q sufficiently small, then, for some 
e > 0, (i) the solution U, and thus C extends to [0,T + e], and (ii) C, is 
hounded and continuous on [0,T + e]. 

Proof. Boundedness and smallness of \U{t)\m on [0, T] follow by Proposition 
14. H provided smallness of C(^) ^-iid iC/oliT"- By Lemma 15.111 this implies 
the existence, boundedness of \U{t)\}j4, on [0, T + e], for some e > 0, and 
thus, by Lemma 14.31 boundedness and continuity of (" on [0, T + e]. D 

Proof of Theorem |j.7| We shall establish: 

Claim. For all t > for which a solution exists with C, uniformly bounded 
by some fixed, sufficiently small constant, there holds 

(260) Q{t)<C2{Eo + ai?). 

From this result, provided Eq < 1/4(71, we have that by continuous induc- 
tion 

(261) C(t) < 2C2E0 

for all t > 0. From (|26ip and the definition of C, in (j259p we then obtain the 
bounds of (|2ip. Thus, it remains only to establish the claim above. 

Proof of Claim. We must show that (|{7| + |C/x|)(^ + V'i+V'2)~"'^ is bounded 
by C{Eo + C(i)^)) for some C > 0, all < s < t, so long as C remains 
sufficiently small. First we need an estimate for U{0,s) and Us{0,s). For 
the inflow case, by boundary condition estimate ()234p and by the hypotheses 
on h{s), we have 

(262) |C/(0, s)\ < C{his) + \U{0, s)|2) < C(i?o(l + s)-' + |t/(0, s)\^) 
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from which by continuity of |[/(0, t)| (Remark 14. 4p and smallness of £"0, we 
obtain a similar estimate to ()261l) : 



(263) 



\U{0,s)\ <CEo{l + s)-\ 



Similarly for an estimate of Ut{0,t), by taking the derivative of (j234p . we 



get 
(264) 



\Us{0,s)\<C{h'{s) + \U\\Usms)) 

<C{Eo{l + s)-^ + \U{0,s)\\Us{0,s)\) 

which by the same argument as above yields 

(265) \UsiO, s)\<C Foil + s)-\ 

Next, for the outflow case with boundary condition (j235p . we have 

\BU{0, s)\ < CEoil + s)"^ + O{\U{0, s)f) 

(266) , 

\iBUUO,s)\< CEoil + sy^ + o{\u\\Usms)). 

Now by (|259p . we have for all t > and some C > that 

(267) \U{x,t)\ + \U^{x,t)\ <C{t){e + i;i+i;2){x,t), 

and therefore 

\Q{U,Uy){y,s)\<CC{tf^{y,s) 
\UiQ{U,Uy)y{y,s)\<CC{tf^{y,s) 



(268) 



with ^ = {e + ^pl + tp2)^ as defined in (f246l) . for < s < t. 

As an estimate for U{x,t), we use the representation (|236p of U{x,t): 



\U{x,t)\ 



+ 



+ 



G{x,t;y)Uo{y)dy 



(Gyix, t - s; 0)BU(0, s) + G{x, t - s; 0)^^7(0, s)) ds 



/ / H{x,t- s;y)UiQ{U,Uy)y{y,s)dyds 
Jo Jo 



Jo 

t roo 



JO 



Gy{x, t - s; y)U2QiU, Uy){y, s) dy ds 



where by applying Lemmas 15. 3115. 61 together with (I268P , we have 
G{x,t;y)g{y)dy 



(269) 



<Eo / {\G{x,t;y)\ + \H{x,t;ymi + \y\ 

Jo 
<CEo{9 + i;i + i^2){x,t) 



-3/2 



dy 
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(270) 



t i-oo 



JQ 



Gy{x, t - s; y)Q{U, Uy){y, s) dy ds 



t fOQ 



<CXitr / \Gy{x,t-s;y)\^iy,s)dyds 
Jo Jo 



(271) 



/ / H{x,t - s;y)UiQ{U,Uy)y{y,s)dyds 
Jo Jo 

<Catf / H{x,t-s-y){e + ^i+i^2?dyds 
Jo Jo 

pt poo 

<Catf / H{x,t-s-y)T{y,s)dyds 
Jo Jo 



and, for the boundary term, we apply the estimate (I263p and Lemma l5.8| 
yielding 

ft 
{Gy{x, t - s; 0)BU{0, s) + G{x, t - s; 0)AC/(0, s)) ds 



(272) 



<C(i^0 + C(i)')(^ + V'l+^2)(x,t) 



for the inflow. Whereas, for the outflow case, noting that G(x, t — s; 0) =0 
in the outflow case, we apply the estimate (12661) . ()267p and Lemma 15.81 to 
give the same estimate as above, yielding 

Gy{x, t - s] 0)5C/(0, s) ds < C{Eo + Citf){e + Vi + V'2)(x, t) 



where we used ([267]) for \U{0,s)\ < C(t)(l + 1)"^/^. 
Therefore, combining the above estimates, we obtain 



(273) 



\U{x,t)\{e + i;i+i;2)'Hx,t)<G{Eo + C{tf). 



To derive the same estimate for \Uxix, t)\, we first obtain by using Propo- 
sition 14.11 



Jo L 

<c{Eo+atf)t-y^, 



dr 



where Bh is the boundary function defined in Proposition 14. IL and thus by 
the one dimensional Sobolev embedding: |C/(t)|yi/3,oo < C|t/(t)|/^4, 

\Q{U,U,U < C{C\t) + AC'El)T 
\Q{U, [/,),,! < C{C\t) + AC^Ei)T 

where T = t-^/^{e + Vi + V'2). 
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Now again applying Lemmas [531l5.8l together with (127411 . (I265|] . and ()266|) . 
we have obtained the desired estimate, that is, bounded by {C'^{t) + CEo){9+ 
^1 + '4'2){x,t), for most terms in the formulation (1242^ of Ux{x,t), except 
one boundary term: 

/ H{x, t - s; 0)\UiQ{U, Uy)y{0, s)\ ds, 
Jo 

which is bounded by C£'o(V'i + 4'2){x,t) by using (j233p . (|267p . and Lemma 

15.71 ^-iid noting that 

\^iQ{U, Uy)y{0, s)\ < C{t)\h{s)\{0 + ^p,+ V'2)(0, s) < CC{t)\h{s)\. 

Therefore, together with (I273p . we have obtained 

(275) {\U{x,t)\ + \U.,{x,tme + %bi+^2r\x,t)<C{Eo + atf) 

as claimed, which completes the proof of Theorem 1 1.7[ D 
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